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We present results for the ground states of a system of spin-1 bosons in a rotating trap. We focus 
on the dilute, weakly interacting regime, and restrict the bosons to the quantum states in the 
lowest Landau level (LLL) in the plane (disc), sphere or torus geometries. We map out parts of 
the zero temperature phase diagram, using both exact quantum ground states and LLL mean field 
configurations. For the case of a spin-independent interaction we present exact quantum ground 
states at angular momentum L < N. For general values of the interaction parameters, we present 
mean field studies of general ground states at slow rotation and of lattices of vortices and skyrmions 
at higher rotation rates. Finally, we discuss quantum Hall liquid states at ultra-high rotation. 

PACS numbers: 03.75.Mn,05.30. Jp,73.43.Cd 



I. INTRODUCTION 

In recent years, there has been considerable progress in 
the art of manipulating cold atoms. By varying experi- 
mental conditions, a number of quantum states of atomic 
matter have been realized. In these developments, an ex- 
citing theme is the parallel between these states of atomic 
matter and states of electronic or vortex matter that have 
been studied before. Two important recent advances are 
the study of spin-full bosons in optical traps and the anal- 
ysis of bosons in rotating traps. 

Optical traps liberate the (hyperfine) spin degree of 
freedom of spin- full atoms, as there is no polarizing mag- 
netic field. This allows a variety of new phenomena, such 
as skyrmions, monopoles and 7r-disclinations. There are 
two different regimes, depending on the sign of the spin- 
dependent interaction c 2 . If this interaction is repulsive 
(c 2 > 0), the system tends to minimize the total spin and 
we speak of the anti-ferromagnetic or "polar" regime. If, 
on the other hand, the interaction is attractive (c 2 < 0) 
the system will tend to maximize the total spin. This 
is the ferromagnetic regime. Examples of such systems 
are spin-1 Bose-Einstein condensates (BEC) which can 
be realized by trapping atoms such as 87 Rb (c 2 < 0) 1,2 
and 23 Na (c 2 > 0) 3 . In both cases, the ratio 7 = c 2 /c 
of the spin-dependent (c 2 ) and the spin-independent (co) 
parts of the (contact) interaction is small, typically a few 
percent, I7I « 0.01 - 0.05. 

Rotating the bosons within the trap leads to the for- 
mation of quantum ground states with a certain amount 
of vorticity stored in the system. Upon increasing the ro- 
tation rate, a single-component (scalar) condensate with 
a repulsive interaction goes through the following stages: 
(i) the nucleation of a single vortex, (ii) the formation 
of a triangular (Abrikosov) lattice of vortices. Theoret- 
ical analysis predicts that this will be followed by (iii), 
a quantum melting of the lattice and the formation of a 
(series of) quantum liquids, where the vorticity is spread 
uniformly over the system. In the presence of internal 



degrees of freedom, such as those associated with the 
spin-states of spin-1 atoms, a similar sequence of quan- 
tum ground states is expected, with additional structure 
provided by the internal (spin) degrees of freedom, and 
by the presence of additional parameters such as 7. 

In this paper we study the ground states of spin-1 
atoms in a rotating harmonic trap, focusing on situations 
where a truncation to the lowest Landau level (LLL) can 
be justified. We consider generic values of 7 in the re- 
pulsive regime cq > 0, and also focus on the special case 
7 = 0, where the interaction has an SU(3) symmetry 4 . 
Before we come to the details of our analysis, we briefly 
summarize some of the results in the literature for rotat- 
ing bosons, with either a single component (scalar case) 
or several components, such as spin 1 (vector case). Much 
of this work is based on the restriction to the LLL. 



A. Non- rotating condensates and topological 
excitations 

The case of vector BEC without rotation has been in- 
vestigated, most often by mean-field theory using the 
spin-1 Gross-Pitaevskii equations, or with the further ap- 
proximation of neglecting the kinetic energy term (the 
Thomas- Fermi approximation). There are two regimes 5 , 
as mentioned already. In the ferromagnetic regime c 2 < 
0, the ground state has maximum possible spin, S = N. 
Such a spin state can be constructed by condensing all 
the bosons into a single-particle state with S z = +1, or 
as a spin-rotation of this. In the opposite antiferromag- 
netic regime, c 2 > 0, the ground state has minimal spin. 
Such a spin state can be constructed by condensing all 
the bosons in the same S z = single-particle spin state, 
or by taking a spin-rotation of this (notice that this is dis- 
tinct from all of the ferromagnetic states). This "polar" 
state breaks spin-rotation symmetry, even though the ex- 
pectation value of the total spin, or of the spin density, 
is very close to zero. The distinct ordered states, that 
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can be mapped onto each other by the (broken) symme- 
tries of spin rotation and phase rotation, are labelled by 
points in an order parameter manifold (or target space). 
For the ferromagnetic case, this manifold is SO(3)^, while 
for the antiferromagnetic or polar case it is S 1 x S 2 /Z<^&. 
These ordered states possess excitations that can be de- 
scribed as topological defects in the order, either with a 
singularity at a point surrounded by a "core" , or without 
a singularity at all. We describe these in more detail in 
the Appendix, but those that carry non-zero vorticity are 
relevant to the rotating case which we discuss next. 



B. Slow rotation: vortices and skyrmions 

If the number of bosons in a rotating trap is sufficiently 
large, the effect of slow rotation can be studied in a mean 
field framework. For a single species of bosons with repul- 
sive interactions, the rotation is accommodated through 
the creation of singular vortices, with vanishing particle 
density at the vortex cores. As the rotation rate is in- 
creased from zero, there is a critical frequency at which 
a single vortex first appears in the system, followed by 
additional vortices at still higher rotation. With a spin 
degree of freedom, the system has several components in 
which to store the angular momentum. There is the pos- 
sibility that a vortex core for one spin component is filled 
by another spin component, leading to core-less vortices 
or "skyrmions". In such configurations, the total parti- 
cle density is nowhere zero, and there is a smooth spin 
texture. Mean field states of this type for rotating spin- 
1 bosons have been investigated theoretically by solving 
the spin-full Gross-Pitaevskii (GP) equations£i2*2i±&. For 
attractive interactions, in contrast, the BEC remains in 
a compact blob, without any vortices, all the way up to 
the maximum rotation frequency (the trap frequency). 



C. Lattices of vortices and skyrmions 

When several vortices are present in a rotating scalar 
boson condensate with repulsive interactions, they line 
up in a triangular (Abrikosov) vortex lattic e) 11 ! 12 . In a 
vector BEC one expects to find similar lattices, built from 
the coreless vortices just described. The details of all this 
depend crucially on the relative strength 7 = c^l cq of the 
spin-dependent interaction. For 7 = 0, where the SU(3) 
symmetry between the different spin components is not 
broken, the lattice that is expected upon rotation is com- 
posed of three intertwined triangular lattices. The vortex 
cores do not overlap, so that the density is (almost) uni- 
form. This lattice has been shown to be independent of 
the strength of the interaction by Kita et al 1 - . The vor- 
tex lattice shows a rich phase diagram, however, when 
the interaction is spin-dependent. For a range of positive 
values of 7, a square lattice composed of 7r-disclinations 
has been predicted^. 



D. Scalar boson quantum liquids 

As the rotation increases, quantum fluctuations be- 
come more and more important and beyond a critical 
rotation rate the vortex lattice is expected to melt. The 
resulting state of matter is disordered, and it has a large 
amount of vorticity stored in it. In this regime, a num- 
ber of quantum liquid states have been proposed. Based 
on the analogy with the physics of electrons in a strong 
perpendicular magnetic field in two dimensions (2D), the 
boson quantum liquid states can be characterized as frac- 
tional quantum Hall liquids. 

An extensive study of the transition to this regime was 
conducted by Cooper, Wilkin, and Gunni^, who per- 
formed exact diagonalization studies in the LLL in a 
toroidal geometry. They predicted that the vortex lat- 
tice melts into a quantum disordered phase at a critical 
value v c ~ 6-10 of the filling factor v — N/N v , with 
N v the number of vortices and N the number of bosons. 
They also found that the quantum state is incompressible 
at v = k/2 < v Cl with k = 0, 1, 2, . . . , and observed that 
the ground states at v = k/2 have substantial overlaps 
with the Read-Rezayi (RR) quantum Hall states^. 

The RR states are incompressible quantum Hall fluids. 
The special case k = 1 is the Laughlin stated for bosons 
at v = 1/2, while the k — 2 state is the Moore- Read 
"Pfaffian" stateA2>. The RR quantum Hall states possess a 
specific "order-fc" clustering property: there is a so-called 
composite-boson order parameteni&Iii, an operator that 
creates k bosons and 2 vortices, and is the minimal order 
parameter (the one that creates the smallest number of 
bosons) that has long-range ( "off-diagonal" ) order. Such 
order implies that the quasiparticles over these liquids are 
(quantized) vortices in the liquid, carrying fractional vor- 
ticity 1/fc, analogous to a fractional magnetic flux &/&o 
that occurs in certain quantum Hall states in electronic 
systems. In quantum Hall liquids, the Hall conductivity 
implies a fundamental quasiparticle charge q = ±j/<I>/<I>o 
in units of the charge of the electron. In the context 
of neutral bosons in a rotating trap, the same argument 
implies that a fractional particle number (q = ±1/2 for 
the RR states) is present in the quasiparticle, relative to 
the background density. Furthermore, for k > 1 there 
are non-local degrees of freedom associated with these 
quasiparticles, that is, the ground states with more than 
three quasiparticles are degenerate in the limit where all 
the separations go to infinity. As these degrees of free- 
dom are of a nonlocal, topological nature, they do not 
couple to local probes and the degeneracy is protected in 
the large separation limit. For the case k — 2 there is 
an interpretation in terms of a Majorana fermion in each 
vortex coraSl. Further evidence for the appearance of 
the Moore-Read state was recently provided by Regnault 
and Jolicoeus2i, who observed the low-lying two-particle 
branch in numerical simulations, upon adding one flux 
quantum in a spherical geometry. They also found evi- 
dence for other quantum Hall states not in the RR series. 
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E. Outline 



Truncation to the lowest Landau level 



In previous work 4 wc analyzed spin-1 bosons in the 
LLL. We identified attractive and repulsive regimes in 
the Cq-C2 plane, and proposed and analyzed two series 
of clustered quantum Hall states (labeled SU(A)k and 
SO(5)k), analogous to the RR states, for spin-1 bosons 
in a rapidly rotating trap. We identified the exact ground 
state for N spin-1 bosons on the disc with one unit 
of angular momentum per particle, the boson-triplet- 
condensate (BTC). 

In this paper we provide further results on the phase 
diagram for spin-1 bosons in the LLL. Employing the 
SU(3) symmetry, we discuss how in a slowly-rotating sys- 
tem with C2 = the exact quantum ground state evolves 
from the non-rotating one towards the BTC at angular 
momentum L = N, and we compute the ground state 
angular momentum L(u>) as a function of the rotation 
frequency ui. Using LLL-mean field theory, we extend 
the results for slow rotation to ^ 0, and discuss the 
various skyrmion lattices. Furthermore, the two series 
of quantum Hall states are discussed in detail and we 
supplement them with a third series. 

This paper is organized as follows. In section [H] we 
define the model by discussing LLL truncation in a 
disc, sphere or torus geometry, specifying the interaction 
Hamiltonian, and make remarks on the general symmetry 
properties for later use. In section lTTTl we study the phase 
diagram by direct numerical diagonalization. In section 
IIVI exact quantum ground state wavefunctions and ener- 
gies for a slowly rotating (angular momentum L < N) 
system in the C2 = limit are presented. For nonzero C2, 
we use a LLL mean field treatment to study the slowly- 
rotating system (in section 0) and the various skyrmion 
and vortex lattices (in section IVTj l. In section IVlTl we dis- 
cuss the quantum Hall states at ultra-high rotation. In 
an Appendix we discuss the topological classification of 
defects. 



In a rotating frame of reference the Hamiltonian for TV 
trapped, weakly-interacting spin-1 bosons is 



II. LLL MODEL HAMILTONIAN AND ITS 
SYMMETRY 



In this section we describe the truncation of the space 
of single-particle states to those in the LLL, and then ex- 
plain the use of different geometries (sphere, torus) once 
this trunction has been made. Then we give the form of 
the interaction Hamiltonian that will be assumed, and 
some analysis of the symmetries of the model, with par- 
ticular reference to certain limits and different geome- 
tries. 



H = 
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Here to is the frequency of the rotation drive, Li the an- 
gular momentum of the z-th particle and H nrt the inter- 
action Hamiltonian, which we discuss below. We have 
set h and the harmonic oscillator length I = (h/rribUJo) 1 ^ 2 
of the trap (with luq the trap frequency and the bo- 
son mass) equal to one. Modes in the direction of the 
rotation axis are frozen out, leaving us effectively with 
a two-dimensional (2D) system. The energy eigenvalues 
of the single-particle part of the Hamiltonian are then 
E n ,rn = (2n + m + 1)u>q — mio, with n > the Landau 
level index and m > — n the z-component of angular mo- 
mentum, labelling the states within each Landau level. 

We consider the model in which the single-particle 
states are restricted to the lowest (n = 0) Landau level 
(LLL)2&. This is valid when the interactions are suffi- 
ciently weak, as we will explain momentarily. The nor- 
malized LLL wavefunctions are <fi m (z) C a with the orbital 
part 4> m (z) tx z" l e~' z ' / 2 (z = x + iy), and Q a a three- 
component complex vector representing the spin state; 
here a labels the eigenstates of the z-component of the 
spin S z for each particle, a =|, 0, J,. [Later in the paper 
it will be convenient also to use the basis of Cartesian 
components for spin 1, labelled by /i = x, y, z.] When 
we use second quantization, we will denote the boson cre- 
ation and annihilation operators for these single-particle 
states by b\ llQ1 b 
numbers by n, 

field operator ip a (z) — ^2 m b ma (j) m {z). The single parti 
cle contributions to the Hamiltonian add up to {loq — lo)L, 
with L = L z i the z-component of total angular mo- 
mentum. We will refer to this geometry as the disc in 
view of the form of the fluid states (for repulsive interac- 
tions) which tend to form a disc or "pancake", because 
of the centrifugal force. Note that we must have ui < lu , 
otherwise the system becomes unstable. 

To study the bulk properties of the quantum ground 
states, we will eliminate boundary effects by using in- 
stead two other geometries and taking the limit ix> — > 
loq. In a spherical geometrji 2 ^, the orbital part of the 
LLL single-particle wavefunctions is 4> m (z) oc z m /(l + 
(J^jr) 2 ) 1+JVl '/ 2 , where z represents position on the sphere 
by stereographic projection to the plane, and R is the 
radius of the sphere. The number of orbitals is restricted 
by the vorticity N v penetrating the sphere, < m < N v . 
The N v + 1 single-particle orbitals form a representation 
of orbital angular momentum equal to N v /2, see Refw 
In the limit R — > oo, keeping N v /R 2 , N and z constant, 
the single-particle wave functions on the sphere reduce 
to those for the disc as above. The total angular momen- 
tum on the sphere is characterized by quantum numbers 



mai "ma, and the corresponding occupation 
bmab-ma- Also, we sometimes use the 
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L for the magnitude, and L z for the z-component. In 
terms of L which has eigenvalues L = ^+ rrii as before, 
L z = ^NN V — L. We emphasize that our definition of L 
when used for the sphere does not have its usual mean- 
ing, but is related to the z-component in such a way that 
the N v — > oo limit agrees with the plane. 

The final geometry we use is the torus. Here the single- 
particle wavefunctions take the form 4>(z) oc f(z)e~ y in 
the Landau gauge, with / a quasiperiodic holomorphic 
function. With N v flux quanta, / has N v zeros in the 
unit cell. There are exactly N v independent solutions, 
of the form f(z) — J\^i $i( z ~ z i\ T )> with T describing 
the geometry of the unit cell and Zi the zeros of /. The 
use of ^-functions ensures that 4> is periodic. Many-body 
states can be classified by their Haldane momentum^. 



B. Interaction Hamiltonian 

In a model description, the Hamiltonian describing the 
2-body interactions of a system of N spin-1 bosons is a 
contact interaction, and contains spin-independent (H n ) 
and spin-dependent (H s ) terms, of strengths cq, c 2 re- 
spectively 

Hi nt — H n + H s 



N 



2^^5( 2 )(r 2 -r J ) 



C + C2 Si ■ Sj 



(2) 



i<j 



Here c = (g +2g 2 )/3, c 2 = (fife— 3o)/3, 9s = ^nti 2 a s /m b 
and as (S — 0, 2) the 2D s-wave scattering phase shift in 
the spin-5 channe&24. A factor 2ir has been extracted for 
later convenience. One can obtain these parameters by 
integrating over the third direction. Assuming, for exam- 
ple, harmonic confinement with quantum length lj_ in the 
^-direction, one finds Ocf = a^ D /V2~jtIj_ when l± <C /. 
For the sphere, the coordinates r in this Hamiltonian take 
values on the surface of the sphere, with radius R. 

The use of the LLL reduced Hamiltonian is justified 
when the interactions are weak. Physical quantities eval- 
uated in the full model differ from those in the LLL model 
by relatively small corrections when the ves <SC 2wq. Here 
v is the typical filling factor (expectation of the occupa- 
tion numbers, summed over a or fj,) of the single-particle 
states. Notice that this condition becomes much less 
stringent as u — > luq in the repulsive regime, as then 
the particles spread out into a pancake, and the filling 
factor v becomes of order 1. 

Finally then, the LLL Hamiltonian in the rotating 
frame which we wish to analyze is 



H u = (cjq - uj)L + Hj 



hit 



(3) 



Note that we use precisely this definition in the case of 
the sphere as well as for the disc. It will be useful also 
to know the ground states of Hi nt for each L. 



C. SU(3) symmetry analysis for C2 = 

In general, the only symmetry in spin space of the 
Hamiltonians H mt and is spin-rotation symmetry 
S , 0(3) sp i n . This implies that spin states will come in mul- 
tiplets of spin S with degeneracy 25 + 1 (with S integer 
since the particles have spin 1). However, at C2 = 0, the 
interaction Hamiltonian reduces to the spin-independent 
interaction H n . In this case the spin-rotation symmetry 
is enlarged from SO(3) sp i n to SU(3) spm . It will be useful 
to understand what this implies about the spin multiplets 
in a finite size system. 

For C2 = 0, the spectrum will contain degenerate spin 
multiplets labelled by S'[/(3)-quantum numbers (p,q). 
These tuples are the Dynkin indices labelling irreducible 
representations of dimension dim( p g j = |(p+l)(g+l)(p+ 
q+2). Since 5*0(3) is embedded in SU(3), each multiplct 
can be decomposed into a set of 5*0(3) multiplets. These 
50(3) spin quantum numbers can be deduced by using 
branching rules for 577(3) 1— > 50(3). The fundamental 
branching rule states that a (p, 0) or (Q,p) multiplet con- 
tains 5 = p, p—2, p — 4, . . . , 1 (0) for p odd (resp., even). 
Using the fusion rule 

(p,0)® (0,«) 

= Me(p-i,5-i)e---e(p-5,o) , (4) 

which is valid for p > q, general branching rules can 
be derived. A multiplet (p, 5) with q odd and p > q 
decomposes in 50(3) multiplets with highest weights 5 
according to the branching rule 



~2~ p+q-2i 

M - s 

i=0 S=2i+1 



(5) 



For q even we find 
(P.?) ~ S )®[ 2 A P° dd ( 6 ) 



9—2 

~2~ p+q—2i 



i=0 S=2i+1 



p+q-2t 



(p,q) h 00 <?je(02ij, p even. (7) 

\ i=0 5=2i+2 ' ^ j=0 ' 

Note that the highest 50(3)-spin in an SU (3)-multiplct 
(j>, q) is always 5 = p + q, and the lowest 5 = or 1. 



D. Orbital symmetry in spherical geometry 

In the plane geometry, i?i nt is invariant under transla- 
tions and rotations in the plane. When working on the 
sphere, this symmetry group is replaced by the rotation 
group 50(3) m b (strictly, we should say 5?7(2) or b when- 
ever N v is odd) of the sphere. In the limit R — > 00 de- 
scribed above, this symmetry becomes translations and 
rotations of the plane. When taking this limit, we also 



D Orbital symmetry in spherical geometry 



5 



C2 






I A 


II 


On = 

^ 






SE/(3) 










Co 


IB 




52 = 



FIG. 1: Overview of C0-C2 plane, with special regions and 
directions marked. 

hold L fixed, and hence many-particle states of definite 
SO(3) oi h quantum numbers (L,L Z ) become in the limit 
infinite-dimensional multiplets of the Euclidean group of 
the plane. States within each multiplet differ only in the 
state of the center of mass variable (which has coordi- 
nate z c — Ylj Zj/N). Thus, if tpL is an eigenfunction of 
Hi n t at certain angular momentum L. then there exists 
a whole "tower" of states ipL+L' oc z^ ip^ with the same 
interaction energy at angular momentum L + L' . 

We remark that in situations where only a few quan- 
tum orbitals are available to the bosons, the spectrum is 
largely determined by symmetry considerations. Partic- 
ular examples are the spectrum for N v = 2 on the sphere, 
where the exact iV-body energies are given in terms of 
Casimir invariants of the orbital and spin-symmetries 
(see eq. (|10|l below), and the case with N v = 4 on the 
torus, where the topological degeneracy eq. $25$ below 
pertaining to particular quantum liquid states is recov- 
ered from the SU (3) spin-symmetry. 

III. MAIN FEATURES OF THE PHASE 
DIAGRAM 

In this section we make a first pass through the phase 
diagram with numerical results on moderate sizes. First 
we consider the ground states of H lrL t in the disc geometry 
for each L, then use this to find the ground states of H u 
as a function of u>. All this has to be done for general 
values of cq, C2. 



A. Global structure of the phase diagram 

First we point out that the magnitude (cq + c|) 1/ ' 2 only 
sets the overall energy scale, so it can be divided out. 
Thus the phase diagram can be thought of as a circle, in 



which a point on the circle represents a ray in the C0-C2 
plane. We wish to examine this for each L, or later for 
each lu. In Figure^the C0-C2 plane is shown with certain 
special directions (co = 0, C2 = 0, go — 0, <?2 = 0) that 
will be important later picked out. 

For L — 0, the ground state has total spin S = N 
for C2 < (ferro regime) and S = (1) for N even 
(odd) for C2 > (anti- ferro regime). These states are 
the way that the broken symmetry states described in 
Sec. 0] (the ferromagnetic and polar states respectively) 
appear in a finite size study. For ci — 0, there is a 
single SU(3) multiplet of spin states, decomposing into 
one SO(3) multiplet of each spin S = N, N — 2, . . . 
The transition at C2 = can thus be viewed as levels 
crossing, with a larger degeneracy on the line C2 = 0. As 
L increases, these two phases at C2 7^ survive in part 
of the phase diagram, as compact drops of fluid, with 
the center of mass carrying all the angular momentum. 
Meanwhile, the positive en axis gradually opens into a 
region that contains other phases. By the time L is > N, 
the C0-C2 plane contains the three regions labelled I- A, 
I-B, and II in figure ^ 

The ground states in regions I-A and I-B are similar 
to the one in the "attractive" regime in the scalar case^. 
The orbital part of the ground state wavefunction is of 
the form Sfr(zi) oc z^. In region I-A (co < 0, C2 > 0), the 
spin state is the same spin-singlet as for the L = ground 
state, and the energy^ becomes [cqN(N — l)/2 — NC2]. 
In region I-B (02 < 0, Co < — c%), the spin state is ferro, 
S = N, giving energy^ (c + c 2 )A(A-l)/2 < 0. At c 2 = 
0, Co < 0, the spin states again form the SU(3) multiplet. 
In the remaining "repulsive" region II, the ground state is 
in general not a common eigenstate of the Co and C2 parts 
of the interaction, and the ground state energy depends 
non-linearly on the ratio 7 = cijc§. Note that we have 
now located the repulsive region more precisely than in 
our previous characterization of it simply as Co > 0. Most 
of the following analysis focuses on region II only, which 
can be parametrized by 7 — C2/C0 alone. 



B. Finite size results in region II as a function of uj 

In figure|21we show the ground state quantum numbers 
(L, S) in region II for N — 6 bosons as a function of the 
rotation frequency lu. As the phase diagram for each lu is 
a circle (which in region II can be parametrized by 7 = 
c 2/ Crjj or by <f> = arctan7), we are free to plot lu radially. 
The parameters are shown in units of luq and with Co = 
0.25, but notice that the ground state quantum numbers 
can only depend on the dimensionless ratios of energies 
(luq — u))/co and C2/C0, so that the structure shown is 
actually present (though with the radial variable rescaled 
and shifted) for all parameter values (unless Co is too 
large). The dashed rays are the lines C2 = — en, c 2 = 
and C2 = co/2, and the outer dashed circle is the locus 
of lu = luq. The ground state angular momentum L and 
spins S at C2 = and as a function of w, are shown in the 
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FIG. 2: Ground state quantum numbers (L, S) in region II 
for TV = 6 spin-1 particles in the planar (disc) geometry, as 
a function of the driving frequency uj (plotted radially) and 
the ratio 7 (corresponding to the angle with respect to the 
horizontal axis). The special directions (72 = 0, C2 = 0, go = 0, 
Co = are shown as double-dotted radial lines. The inset 
shows a cut along the C2 = direction, with the angular 
momentum given on the vertical axis and the (degenerate) 
spin values S marked at each of the steps. In this figure, the 
parameters Co, C2, and u) are in units of Wo, and the value 
Co = 0.25 is used in the main figure as well as in the inset. 
For additional discussion, see the main text. 



inset. The degenerate spin values at L < N are seen to 
correspond to the following irreducible SU(3) multiplcts: 
(p, q) = (6, 0) for L = 0, (p, q) = (4, 1) for L = 1, (p, q) = 
(2,2) for L = 2, (p,q) = (1,1) for L = 3 and (p,q) = 
(0, 0) for L = 6. Note also that for C2 < 0, Co > — C2 
the ground state spin gradually decreases from S = 6 at 
w = to 5 = 0. 



IV. SLOW ROTATION: EXACT GROUND 
STATES AT c 2 = 0, L < N 

For larger sizes, a brute-force numerical approach is 
not feasible, so we develop other approaches. In this sec- 
tion we determine the exact ground state energies and 
wavefunctions for slow rotation (angular momentum up 
to the boson number, L < TV), for c > and C2 = 0, ex- 
ploiting the SU(3) symmetry described in the section ITU 
Some of the ground states we find were described in 
Refill. We analyze a system of N spin-1 bosons in spher- 
ical geometry with N v quanta of vorticity, with the disc 
geometry emerging as the limit N v — > 00. We remark 
that for N sufficiently large, it becomes natural to dis- 
cuss low energy properties in terms of mean field con- 
figurations that break the various symmetries and whose 
energy is slightly higher than that of the exact quantum 
ground state; this will be discussed in section Ivl 



A. Exact eigenstates of H n 

The ground state spectrum for C2 = and L < N can 
be understood by exploiting the SU(3) symmetry of the 
Hamiltonian H n . In our analysis we proceed as follows. 
We consider two series of eigenstates of H ni in which 
(roughly speaking) the bosons occupy at most the lowest 
three orbitals. Among these eigenstates, we identify the 
exact quantum ground states on the disc and the sphere, 
as a function of the angular momentum. This then al- 
lows us to compute the to dependence of the ground state 
angular momentum for general N at C2 = 0. 

We write the first series of eigenstates as \p, q, n) 1 . 
These states contain doublets and triplets of spin-1 
bosons that are fully antisymmetric in spin indices, and 
in the orbital indices (guaranteeing the overall symmetry 
that is required). The different numbers of single bosons, 
doublets and triplets correspond uniquely to the values of 
N and the quantum numbers (p, q) of the corresponding 
SU(3) multiplcts. The triplets, which appear n times, 
are singlets under SU(3), and so do not affect the over- 
all SU (3) representation. The highest spin component 
(S z = p + q) of the corresponding SU(3) multiplet takes 
the following form (up to normalization) 

\p,q,nY cx 

[el • B\f[e 2 ■ (B\ x ■ (§\ x bJ)]«|0), (8) 



with ei = (1,0,0), e 2 = (0,0,1) and B\ = 
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P , Oo''i,«i''L)' Clearly, the total number of bosons is 
N = p + 2q + 3n. The energies corresponding to eq. © 
are 

^p.g.n/co = a? v n(n-l) + ot$ v q(q-l) + ~p(p-i) 



N " N 

+a 3 v np + -gp -I- a 4 v nq, 



(9) 



with 



_ o HJV^-20iV„+6 
~~ "*4(2iV„-3)(27V„-l) 

v A f „ _ 7jV„-4 
l 3 — 2(2JV„-1) 



N v _ 5Ny-2 
'2 ~ 2(2JV„-1) 



, jV„ 5a, " 



a:. 



This energy is for spherical geometry and it depends on 
the number N v of flux quanta. For N v —> oo eq. 
gives the energy in a disc geometry; N v — 2 gives the 
energy on a sphere with 3 orbitals. On the basis of exact 
diagonalization studies for N = 6,9,12,15,18 particles 
we claim that on the disc for L < N/2, the ground state 
multiplet is precisely \p, q, 0) 7 , with p — N — 2L, q = L. 

On the sphere with N v = 2, we have obtained a much 
stronger result, namely a closed form result for all 
eigenvalues of H n . It turns out that these energies can 
be given in terms of the number N of bosons, the total 
angular momentum L and the (p, q) labels of the SU (3) 
representation, according to 



t?N v =2 I 



CO 



18 



N(N — 1) + 



-T 2 

6 p ' q 



-L(L 
6 v 



1) 



(10) 



where T 2 q = 



(p 2 + q 2 + pq)/3 + p + q is the quadratic 
Casimir operator for SU (3) in the representation (p, q) . 
Specializing this expression to the states in scries I, by 
eliminating N in favor of n and using the fact that L = 
p + q, reproduces the result in eq. J§| for N v = 2. 

Analyzing the ground state on the disc for L > N/2, 
we identified a second series of states \p, q, n) H . One can 
think of the type II states as having the p single bosons 
in to = 1 rather then m = 0, so that now e\ — (0, 1,0). 
That is not quite correct for the energy eigenstates, as we 
will explain below, but it does give the correct quantum 
numbers. The states in series I, II share the property of 
having p single bosons and q doublets, leading to SU (3) 
Dynkin labels (p, q) . It may be illuminating to display 
the structure of the states in terms of diagrams similar to 
Young tableaux as in figure |3| For the orbital structure 
of the highest-weight states in either series I or II, the 
lengths of the three rows represent the number of bosons 
in the orbitals m = 0, 1, 2 respectively (in the rough 
point of view, which will be corrected below), while the 
differences p, q and n in the lengths correspond to the 
SU(3) structure. Essentially these diagrams are ordi- 
nary Young tableaux for the states, but with the first 
two rows exchanged in the case of series II. 

For the case of the type II states, the following cor- 
rection must be made to obtain the energy eigenstates. 
In the case of scalar bosons, it is known^LSJiSi that the 



•a 



D-D 



\p,q,n}' 



\p,q,ri)' 



FIG. 3: The structure of the two different series of eigenstates 
of H n , displayed in a form similar to Young tableaux. In both 
cases, the corresponding SU (3)-representation has Dynkin la- 
bels (p,q). 



ground state configuration at L = p of p bosons is a vor- 
tex located at their center of mass, with wavefunction 
Jl^Zi — z c ) with z c = J2i Z i/P- This state is not entirely 
restricted to the m = 1 orbital, as there are components 
in which other orbitals in the range < m < p are oc- 
cupied as well. The p bosons in the state \p, q, n) n form 
such a vortex. This complication makes it difficult to 
write down the closed form expression for the states in 
series II; based on numerical analysis for small N and 
mean field results for large ./V (see section 0, we do pro- 
pose the following closed form expression for the corre- 
sponding energy on the disc 



E p.q.n/ c 



33 5 1 

—n(n - 1) + -q(q - 1) + -p{p - 2) 



25 11 

+ ~7T n< l + ~7T n P + IP- 



(11) 



Note that the p-independent terms in this formula are 
identical to those for type I states with N v = oo. The 
state \p, 0,0) 77 , has energy p(p — 2)/4, which is exactly 
the ground state energy of a rotating scalar BEC at L = 
p = N. This justifies the interpretation of the polarized 
subsystem with p bosons forming a vortex at the center 
of mass. However, it turns out that \p, 0, 0) 77 will never 
be the lowest energy configuration for a rotating spin-1 
system. 

Among the type I/II states the following are special. 
First, \p, 0,0) 7 is the non-rotating ground state, corre- 
sponding to the (p, 0)-multiplet. Second, |0, q,0) gives 
a wavefunction composed of anti-symmetrized pairs of 
bosons, a Boson-Doublet-Condensate (BDC) or (0, q)- 
multiplet. Third, |0,0,n) is composed of 3-body singlets. 
It is a condensate of triplets or boson-triplet-condensate 4 
(BTC); we shall see that it forms the ground state at 
L = 3n = N. The BTC-state can be regarded as a 
symmetrized version of the core-less vortices observed in 
mean field studies (see section [V] for more on this). 

More generally, the type I/II states are examples of 
"(multi-) fragmented" condensates^, see also2&, in the 
sense that they contain several macroscopically occupied 
elements in the density matrix. For instance, for the 
BTC- and for (any component of) the BDC-state we 
have (77 mQ ) B DC = (1 - <W)(1 - S al )q/2, (n ma ) BTC = 
n/3. Since the spin is fixed in these states, (At7 q ) 2 = 
((t7 q - (t7 q )) 2 ) = 0, where 77 Q = J2 m =0,i,2 n ma- However, 
within each spin component, the fluctuations of the boson 
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UJ 

FIG. 4: L/N as iV — *• oo of the ground state on the disc 
(bold line) and on the sphere with N v = 2 (dashed line), as 
functions of <Z> = (u> — ujo) / (coN) at C2 = 0. The horizontal 
lines mark the values L/N — 1/3 and L/N = 2/3. The cusps 
in both curves indicate the point where the m — 2 quantum 
orbital is first used in the ground state. 



number between orbitals is of the order of the system size: 
(An mQ )| D c = q(q + 2)/12, (An mQ )| TC = n(n + 3)/18. 
This is an indication that, as in the case of the singlet 
ground state at L — in the antiferromagnetic regime 
and the related "polar" mean field stated, it may be 
best to think of these states as broken symmetry states^. 
That is the approach we will take in Sec. 



B. Exact ground states at C2 = as a function of L 

or ui 

The ground state of a rotating gas with N spin-1 
bosons in the LLL and a spin-independent (ci = 0) in- 
teraction is formed by a sequence of type I or II states 
lying on a certain path in (p, <7)-space as L increases. To 
find the ground state in a rotating frame of reference, 
we need to find the ground state of H w , eq. ©, instead. 
Since this Hamiltonian contains only two energy scales, 
the ground state angular momentum per particle L/N 
can be written as a function of the ratio (uj — ujq) / (N cq) . 
For finite boson number this function consists of a se- 
quence of steps, as can be seen in figure |21 It turns out 
that (thanks to our judicious choices of factors of N) the 
limit N — ► oo with L/N and (uj — uio)/(coN) fixed of this 
function exists, and this is the most convenient informa- 
tion to display. In the following we determine the path of 
the ground states in (p, q)-space as a function of L, and 
the L(uj)/N behavior of the ground states in this limit 
for both the sphere (N v — 2) and the disc (N v = oo) in 
the regime L/N < 1. 



1. Ground states on the sphere at N v = 2 

On the sphere, our notion of rotation is such that the 
50(3)orb-angular momentum L decreases as the system 
rotates faster and faster. With three orbitals (N v = 2) 
we have L z = N — L (see section III All . (We consider 
N v = 2 because this case can just accommodate L < 
N.) At L = N, we know already that the |iV, 0, 0) 7 
multiplet forms the ground state. As L starts to decrease, 
again a type I state has the lowest energy; the (p,q)- 
path is parametrized by (2L — N,N — L). Bosons are 
gradually added to the m = 1 orbital and form anti- 
symmetrized pairs with the remaining ones. The point 
up to which this continues can be found by comparing 
the energies of \2L - N, N - L, 0) 1 and \2L — N + t,N — 
L — t, t/3} 1 . After minimizing with respect to t this yields 
the critical S"[/(3)-indices (p,q) c = (N/3, N/3). At this 
point, with L = 2N/3, ground states with a nonzero (n > 
0) number of triplets become energetically favorable. In 
the remaining region, 2iV/3 > L > 0, type I states are 
the ground states following the path (p,q) = (L/2,L/2). 
Eventually this terminates on the BTC at L = 0. L/N 
of the ground state as a function of the rotation drive u> 
shows a cusp at L/N = 2/3 (L/N = 1/3), as is shown in 
figure |U 



2. Ground states on the disc 

For a system on the disc (N v — oo), the results are 
rather different. We will again present the ground states 
in order of increasing L. At L = 0, the \N, 0, 0) J - 
multiplet forms the ground state as we know. For 
L < N/2 the ground state is formed by a type I state with 
n = and SU (3)-quantum numbers (p, q) — (N — 2L, L). 
This state terminates on the BDC at L — N/2. In 
this range, increasing L leads, as on the sphere, to 
more bosons occupying the m = 1 orbital, forming anti- 
symmetrized pairs with the ones in the m = orbital. 
For L > N/2, the type II states have the lowest energy. 
As L increases, bosons move from the m = into the 
m = 1 orbital, decreasing the number of doublets, and 
giving type II states at (p, q) — (2L — N, N—L). Compar- 
ing the energies of \tN, (l-t)N/2,0) n and \(t-s)N, (1- 
t)N/2, sN/3} 11 , we can determine the point where it be- 
comes favourable for triplets to enter the ground state. 
We find a critical angular momentum L — (1 — t c )N with 
t c r~j 1/3 — 3/N, which approaches L = 2N/3 for N large. 
For L > 2N/3 the number of triplets is gradually increas- 
ing as L grows. Minimizing \2L — N — s, N — L, s/3) 11 
with respect to s, we find that the ground state is 
now the type II state with s(N,L) = 3L — 2N, giving 
(p, q,n) = (N — L,N — L,L — 2N/3). For L ~ N the 
ground state is the BTC with p = q = 0, n = N/3. 

To summarize the above results, for < L < N/2 the 
ground state is given by \N — 2L,L,Q) 1 and for N/2 < 
L < 2N/3 by \2L-N, N-L, 0) n . In the remaining range 
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2 A/3 < L < N the number of 3-body singlets is nonzero, 
and the ground state is given by \N — L, N — L, L — 
2N/3) 11 . Minimizing the energy in a rotating frame of 
reference leads to the L(w)/A-dependence of the ground 
states for A — ► oo which is depicted in figure 0] In this 
figure, the curve shows a cusp at the point where the 
to = 2 orbital first enters the ground state configuration, 
which is at L/N = 2/3 for the disc with A large. A 
signature of this cusp in an experimental system might 
be a change in the expansion rate (the rate of change 
of the outer radius of the drop with respect to lo) if the 
angular momentum exceeds 2A/3. We shall see that the 
cusp survives in the anti-ferromagnetic regime, C2 > 0. 

It is important to contrast all this with the well-known 
behavior of scalar bosons in a rotating trapSiSS. In the 
latter case there is a jump from L/N = to L/N = 1 (for 
all N) when one vortex enters the system, whereas for 
spin-1 bosons we find (at N — > oo) a continuous L(u)/N 
curve with a discontinuous slope. 



V. SLOW ROTATION: LLL MEAN FIELD 
THEORY 

At low rotation rates, the typical boson occupation 
numbers (n ma ) of the occupied (n ma ^ 0) single-particle 
states are large compared with 1. In this situation, a 
mean field (or classical) approach to the problem is gen- 
erally expected to be quantitatively accurate. In such an 
approach, the boson operators are replaced by expecta- 
tion values, which are complex c-numbers: b\ na — > b* mal 
and the second-quantized Hamiltonian is then minimized 
with respect to both the magnitude and phase of these 
numbers to find the ground states. In essence the re- 
sulting state is a Bose condensate with the bosons con- 
densed in one linear combination of the single-particle 
states. This typically involves breaking the orbital and 
spin symmetries, as well as particle number conserva- 
tion. (States with definite values of the good quantum 
numbers such as N, S, L can be obtained afterwards by 
applying a projection to the mean field quantum stated.) 
In the case of very low rotation, where L < N, we have 
seen that (neglecting the subtleties that arose for type 
II states) the states essentially involve only the to — 0, 
1, and 2 states, so that the basis set for the mean field 
calculation is particularly small. In these cases, the mean 
occupation numbers of the single particle states are of or- 
der N, and their energies exceed the exact ground state 
energy (which is of order TV 2 ) by an amount of order 
N. In this section we pursue this mean field calculation 
for this regime. This gives us easy access to the ground 
states at large N for C2 ^ in region II. In the following 
section, we study instead the mean field states at larger 
rotation, which can be assumed to be states in which the 
translational and rotational symmetry group of the plane 
is broken to that of a lattice. 

In terms of the complex numbers b ma , b^ a , the energy 
becomes a quartic polynomial and the ground state can 



be found by minimizing this polynomial with respect to 
these variables. This is done here with the mean boson 
number (N) = J2 m a Kncfima and angular momentum 
(L) = Ylma m ^mofima fixed at the values A, L, respec- 
tively. The spin is not constrained at all. For C2 = the 
Hamiltonian on the sphere takes the form 

H n = Cp ^ ^ ^ t Vm 1 m 2 m3m i b mia b rn2 pb m3a b m4 f3, 



ct(3 mi-nu 

with matrix elements 



V m 



(12) 




2N V 
m 3 + to 4 



rmi+m 2 
m 3 +m 4 • 



(13) 

This exhibits the dependence on the spatial orbitals. For 
c 2 =/= 0, the matrix elements in the additional term consist 
of V mi ... mi multiplied by matrix elements of Sj • Sj , which 
depend on the ot\, . . . , 0:4 labels of the bosons. These 
matrix elements can be found in standard quantum- 
mechanics texts. 

The fact that mean-field configurations break the vari- 
ous symmetries implies that the minima of the mean field 
energy form orbits under the action of these same sym- 
metries. On a disc, and at C2 ^ 0, one expects and finds 
that, typically, from a generic minimum there are 5 flat 
directions leading to adjacent minima with equal energy. 
These flat directions correspond to the 3 generators of 
the SO(3) spin symmetry, an overall phase, and an or- 
bital 0(2) rotation. For spin-independent interactions 
the symmetry orbits are generically 10-dimensional. 

One convenient quantity to plot is the expectation 
value (S) of the spin, whose length is conserved under 
global spin rotations. In special cases, this expectation 
value is axi-symmetric; in the more general case it is non- 
axisymmctric and the mean field configuration breaks the 
orbital 0(2) symmetry. Another useful quantity is the 
three-component condensate wavefunction (analogous to 
the familiar spinor for spin-1/2), which is the expectation 
value of the field operator, (ip a (z)) = J2 m b m a<l>m(z) (see 
Sec. Ill A|) . It is a vector in the a =f, 0, J. basis. From 
this we can plot the density in each spin component in 
position space. This could be accessed experimentally if 
after switching off the trap to allow the particle cloud to 
expand, a Zeeman term is switched on, which causes the 
three a components to separate as they expand. 

As an example, we plot in figure [5J the 2D density pro- 
file in each spin component of two different mean field 
ground state configurations at C2 = 0, L = N. The top 
frame shows the densities for the condensate proportional 
to (4>o, (f>i,4>2)', the lower frame shows a configuration that 
is related to this by an 5/7(3) rotation. The total density 
in each of the m = 0, 1, 2 orbitals is an SU(3) invariant, 
and it is the same for both configurations shown in fig- 
ure [SJ The mean field energy of these configurations is 
Emf = 3g A 2 , in agreement with order N 2 term in the 
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FIG. 5: Two-dimensional density profile of each of the spin 
components of two LLL mean field ground state configura- 
tions at C2 = 0, L = N. The upper panels shows the 
axi-symmetric spin-vector configuration ((j>o[z),cj>i(z),(j)2(z)). 
The two configurations share the same distribution of the to- 
tal density, and they are related by the SU(3) symmetry. For 
C2 7^ 0, there are similarly distinct profiles related by SO(3) 
symmetry. 



***** 




FIG. 6: Regions in the 7, I plane in which only the m = and 
m = 1 orbitals are present in the mean field ground state on 
the disc are shaded. The angular coordinate is <j> = arctan7 
and I = L/N is plotted radially. In the shaded strip near 
I = 1, a "polar vortex" forms the ground state. The dotted 
lines mark the 7 = ±1 and 7 = directions. 



energy of the exact quantum (BTC) ground state, eq. @ 
with N v = 2, p = q = and n = N/3. 

First we consider the disc geometry with ci = 0. Car- 
rying out the mean field minimizations, we find in terms 
of £ = L/N that for < £ < 2/3 the number densi- 
ties (n m ) = J2 a Kncfima in the orbitals of the mean field 
ground states behave like (here and in the remainder of 
this section, these numbers are normalized so that they 
sum to 1) (no) = 1 — £, (ni) = £ and (n?) — 0. For 
\ < £< 1 we find (n ) = ± (m) = |-£and (n 2 ) = 



All this is in agreement with the results derived from the 
exact quantum ground states in section HVI 

For very small interaction ratios I7I <§; 1, the total 
densities in the orbitals remain the same as for 7 = 0, 
but there is non-trivial structure in the spin dependence, 
leading to spin transitions at critical values of I = L/N, 
as we will describe shortly. 

In figure HO we have plotted region II of the phase dia- 
gram, this time with £ radially. The shaded regions show 
where only the first two orbitals (m = 0, 1) are present in 
the condensate. One region is a tiny strip near 1=1 for 
7 > (7 + 4\/2)/17 w 0.75, where the (m,a) = (1,0) state 
is occupied by all the bosons. This state can be seen as 
a "polar" vortex, since it has the same spin state as the 
polar BEC. The other region, centered (roughly) around 
the C2 = axis, contains states in which both the m = 
and m = 1 orbitals are used. 

In the anti-ferromagnetic regime for £ < 1 there is a 
large area where the m = 3 orbital requires a non-zero 
density; in this area, mean field theory in which only the 
first three orbitals are used is not valid. However, around 
and at the 5J7(3)-axis and around the polar vortex as well 
as in the ferromagnetic regime, the density in the m = 3 
orbital is very small for £ < 1 and can safely be ignored. 
Besides, if the energy H u in a rotating frame (see eq. ©) 
is minimized, only the states which use the first three 
orbitals m = 0, 1, 2 arc of interest for £ < 1. [This is with 
the exception of the vicinities of the boundaries of region 
II (see figure ^1 at 7 — » 00 and at 7 = — 1.] 

In the following subsections we present results for the 
LLL mean field ground state for I7I <§; 1, in both the fer- 
romagnetic and anti-ferromagnetic regimes, and we dis- 
cuss the ground states at I = 1 for general values of 7. 

Our mean field results pertain to the LLL, relevant 
for the regime of weak interactions, and they thus differ 
from the mean field solutions of the GP equations^^^. 
Nevertheless, there is agreement on some of the impor- 
tant features, such as the smooth dependence of L on uj 
in the ferro regime, and the role of the state with a sin- 
gle 7r-disclination near £ = 0.5 in the antiferromagnetic 
regime^*. 



A. Anti-ferromagnetic interactions 

We now specify the mean field ground states, given in 
the form of a three-component condensate wave function, 
for small, positive 7 = +e, and for £ < 1. As before, the 
condensate wave function is a vector in the a =|, 0, J, 
basis. In the table below we specify the mean occupation 
numbers of the four states that we found. 

Note that the condensates given in this table are spe- 
cific representatives of families of condensates that are 
related by the SO(3) spm symmetry. There are two crit- 



ical values, £i 



V2 
7 



0.37 and £% = 10 - 4\/5 - 
W85 - 38\/5 « 0.83, where we see a discontinuous re- 
arrangement of the condensate configuration and of (S). 
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TABLE I: Occupation numbers of the LLL mean field ground 
state with small anti-ferromagnetic interaction 7 = e. 



For 7 = +e, 2/3 < ^ < i\, the condensate can be 
represented as (rtfii (z), afoiz), £<fio(z)), with £ = 

< £ < 1 we 



= t/L- 



2iY 



4JV 
3 



3 > ' — v 3 L, while for c b 
have (-a(j> 2 {z),T(f) 1 (z),£ i (l) Q (z)). 

In figure |H1 we have depicted the ground state angular 
momentum per particle, £, as a function of the rotation 
frequency u for some positive values of 7. It is seen that 
upon increasing 7 a semi-plateau (a distinguished part 
of the curve on which the angular momentum increases 
gradually) develops. Upon increasing 7 further, the semi- 
plateau becomes flatter and the width decreases, until for 
7 larger than some critical value 7 C ss 1.19, £(uj) jumps 
from I = to an t — 1 plateau at a critical frequency u> c 
given by u>q — w c ~ 0.15co-/V. This is a transition from 
the non-rotating state to the polar vortex, analogous to 
what occurs in the scalar boson case. 



FIG. 7: Two-dimensional plot of (S) at both sides of the spin 
transition at £ = £%. The intensity codes the length \{S)\, 
while the color indicates the direction on the spin sphere as 
in figure [TU1 below. The left and right pictures correspond to 
eqs. 1141 and I15L respectively. 



B. Ferromagnetic interactions 

With small negative 7 = — e the mean field ground 
states for slow rotation are characterized (up to 5'0(3) S pi n 
rotations) by the occupation numbers given in table M 
Again, we find two spin transitions, the first at £~ e = 
2 - V2 w 0.59 and the second at £7 e w 0.69. 



For nonzero 7, these changes in the condensate are con- 
tinuous; they become singular (discontinuous) only as 
7->0+. 

For £ < £ e a , the condensate can be represented by 

(^M*)ivM*)>-j3M*)) with A = Vn^l, t] = Vl. 

Applying S , 0(3) sp in rotations, one finds alternative rep- 
resentations such as (-^[\(j) (z) — r/(f)i(z)], 0, -^{\(j) (z) + 

n<f)i(z)]). The S , 0(3) sp i n -invariant quantity | (-S") | 2 is 
found to be 



1(5)1 



JV 2 



£{\-£){z + z) 



z) 2 e- 2 ^ 



The state that emerges at £ 
(i](j) 1 (z),0,X(f>o(z)), leading to 

N 2 



> a 



(14) 

corresponds to 



K^^ai-^M 2 ) 2 ^ 



(15) 



For £ 1/2, the integrated value of (S) for this state is 
non-zero and there is a spontaneous magnetization. In 
figure a two-dimensional plot of (S) at both sides of 
the spin transition at £ = i e a is shown. The state at 
< £ < £ e a can be viewed as a configuration of two ir- 
disclinations off the center of the trap, while the state in 
the regime £ e a < £ < 2/3 (or possibly even as far as £1) 
can be understood as a single 7r-disclination in the polar 
state. 

The angular momentum for which the m = 2 orbital 
is first occupied in the mean field ground state, £ = |, 
is robust against small anti-ferromagnetic interactions. 
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TABLE II: Mean occupation numbers of the LLL condensate 
for small ferromagnetic interaction 7 = — e. 



For £ < £~ e the condensate can be represented by 
(0, rj(j)i(z), A0o(z)) with A and 77 as given above. In this 
state, the expectation values of the components of the 
spin vector take the following form 



(S x 



(Sy 



N 

N 
7rV2 

N 



y/£(l-l)(z + z)e 



-W 2 



y/e(i- e)(-i)(z-z)e 



-\A 2 



(S z ) = -(l-£)e 

7T 



-W 



(16) 



The state at £~ e < £ < 2/3 corresponds to 
(0, \4>o(z), 7](f)i(z)), leading to a spin vector that vanishes 
at the center of the disc. The spin-textures for 7 = — e, 
£ < 2/3, can be interpreted as half-skyrmions (or merons) 
(see also section V D below) . 



For 7 



-e, 2/3 < £ < £ b 



the condensate can 



be represented as (<t4>2(z),^<Po(z),t4>i(z)), with £, a 




FIG. 8: The ground state angular momentum per particle I on the disc as a function of u> = (oj~ uio)/coN for various interaction 
strengths and slow rotation. Upper figures: anti-ferromagnetic regime, <f> = .1, .5, .75. Lower figures: ferromagnetic regime, 
4> = -.l, -.3, -.5. 




FIG. 9: The spin texture at I = 1 for 7 = ±e. The left side 
is the anti-ferromagnetic case, the right is the ferromagnetic 
case. The color coding is as in figure ED 



and t as above, while for £ b £ < £ < 1 we have 

{c7(h{z),T<t>l{z),t,<fo{z)). 

In the ferromagnetic regime the lu dependence of the 
ground state angular momentum becomes a smooth 
curve; see Figure |H1 



C. Mean field configuration at L — N 

Assuming that only the first three m = 0,1,2 or- 
bitals participate in the ground state, we find that 
the mean field ground states at £ = 1 take the form 
(T&4>1 (z) , T<f>i (z ) , £(j> (z) ) , with £ = a = y/k^ and 
t = x/l — 2fc±, and with the (+)— sign corresponding 
to (anti-)ferromagnetic interactions. The parameters k± 



depend on 7 according to 

T(19 + 28V2)7 2 + (42 ± 4\72) 7 - 3 



k± (7) = 



71 7 2 + 1267 - 9 



The orbital occupation numbers, given as 
(n 1Q ) = 1 - 2fc±, (n ;) = (n 2 f) = fc± 



(17) 



(18) 



are continuous for 7 going through 0, but the spin- 
texture, which is sensitive to the phases in the condensate 
wave function, is not. We find that for 7 = ±e, up to an 
overall constant, 



(Sx) 
(Sy) 



= (1 T ^=zz)(z + z)e-^ 2 

= (l T -Lzz)H)(z-z)e-M 2 

= (l-i(2*)V W ". 



(19) 



Note that in the anti-ferromagnetic case, the expecta- 
tion value of the spin vector is vanishing on the circle 
zz = V2, while in the ferromagnetic case we see a sin- 
gle skyrmion texture with (S) non-vanishing everywhere. 
Figure shows the spin texture at £ = 1 for 7 = ±e. 

From (|17fl it is possible to derive the critical anti- 
ferromagnetic interaction ratio for which the polar vortex 
appears, by simply solving £+(7) = 0. The critical value 
found then is 7* = (7 + 4\/2)/17 w 0.75 (see figure©. If 
7 increases towards 7*, the density in the m = 3 orbital 
acquires a small value. So, strictly speaking, the states 
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discussed here are not the true mean field ground state in 
the whole intermediate region. Around 7 = and 7 = 7* 
however, (n^ a ) is zero and the value of 7* is in agreement 
with numerical results. 

In the ferromagnetic regime, upon lowering 7 the pa- 
rameter k- gradually decreases from fc_ = l/3at7 = 
to fc_ = 1 — ^= at 7 = — 1, with the corresponding occu- 
pation numbers given in eq. (|18fl . 



sate are given in table IIHI (up to 50(3) sp i n and SO(3) m b 
rotations). In the trajectory from £ = to £ = 1 
there are no spin transitions. The I — 1 state, which 
has (iioi) — (nio) = (^2|) — is the mean field 
ground state for arbitrary ferromagnetic spin interac- 
tions, > 7 > — 1. It is a single skyrmion texture with 
both uniform number density and magnitude of the spin 
density, and is discussed further in the Appendix. 



D. The sphere with N v = 1, 2 

It is instructive to perform LLL mean field theory on 
a system of spin-1 bosons in a spherical geometry, with 
N v = 1 or N v = 2, meaning that 2 or 3 orbitals are avail- 
able to the particles. To compare with the disc as before, 
we write these results in terms of £ = (| JV„ — L) /N. No- 
tice, however, that by flattening out the sphere by stere- 
ographic projection, the results are qualitatively similar 
to those for the disc when only the first two or three or- 
bitals are occupied. This is especially true for the states 
at N v = 2, £ = 1. Even though the topological classifica- 
tion of textures (see Appendix A) does not strictly apply 
to the plane, the form of the spin textures on the sphere 
is a useful guide to those in the disc for I < 1. 

For the case N v = 1 (two orbitals on the sphere), we 
mention the following results. With < 7 < 7r/4 the 
ground state configuration is the same as the one we 
found on the disc for £ < £ e a . This configuration can 
be interpreted as two 7r-disclinations at opposite poles 
of the sphere. For 7r/4 < 7 < 7r/2 all bosons occupy 
the a = spin component, forming a polar state with 
a single vortex. In the ferromagnetic regime, with very 
small 7 we find the same spin transition as the one on 
the disc at £ — £~ e . With N v = 1 this transition lies at 
(. = 1/2. These configurations can be interpreted as a 
half-skyrmion (or meron) in the spin texture, with the 
spin density vanishing at one point on the sphere, around 
which the spin density winds around the equator in S- 
space, passing over one pole at the opposite end of the 
sphere. If the interaction is deformed by increasing N v 
towards N v — > 00, the location of the spin transition is 
gradually shifted towards £ = £~ e . With finite ferromag- 
netic interaction (N v = 1 again), there is a finite region 
where the core traces a path over the sphere (from the 
south pole to the north pole, as £ increases) and connects 
the two sides of the transition. The interaction energy 
is clearly independent of £. This region is bounded by 
j(£) = -|arctan(2<? - 1)| for < £ < 1. 

For N v — 2 and 7 = 0, the occupation numbers, 
summed over spin, in the three available orbitals are 
given by (n ) = 1 — £, (m) = £ and (n 2 ) = for < £ < 
|, followed by (n ) = f - \£, = |, (n 2 ) = \£ - \ 
for I < i < 1. These mean field results agree with the 
exact quantum ground state results obtained in section 

M 

For N v — 2 and small ferromagnetic interactions, 
7 = — e, the mean occupation numbers in the conden- 
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TABLE III: Mean occupation numbers of the LLL mean field 
ground state in spherical geometry, N v = 2, with small ferro- 
magnetic interaction 7 = — e. 



In the case N v = 2, and small anti-ferromagnetic in- 
teractions, 7 = +e, for < £ < | the mean occupation 
numbers per orbital of the condensate are the same as in 
the ferromagnetic case, but the spin structure is different. 
For £ > i the spin expectation values in the m = and 
2 orbitals become non-zero (without a discontinuity) and 
are not linear functions of £. Since | (S m =i) | 2 = and the 
number density is constant in the m = 1 orbital, the spin 
state describing the bosons in this orbital can be arranged 
by an SO(3) rotation to be (0, 1, 0)/-\/3. The vectors 
representing the bosons in the m — 0, 2 orbitals then 
are simply constructed. Together with the previously- 
mentioned vector they form a mutually orthogonal set 
which minimizes H n . Provided that the spin- vectors 
are properly normalized, the energy can be expressed in 
terms of one parameter a(£), which is connected to the 
spin densities by cos[2a(£)\ = \{S )\/{n ) = \{S 2 )\/(n 2 ). 
Minimizing the energy with respect to a(£) gives 



a{£) 



\ 



1 - 4A + 6A 2 



(20) 



with A = i(£ — -|). The maximum of the anti- 
ferromagnetic energy is not dependent on the angular 
momentum and lies at a = n/2. In the ferromagnetic 
case this point minimizes the energy, corresponding ex- 
actly to the occupation numbers in table ITTT1 

At £ = 1 , there are solutions with uniform density, with 
an unbroken 50(3) subgroup of the S , 0(3) or b x 5 , 0(3) sp i n 
symmetry, as the limiting case of the previous £ < 1 
states. This case is also discussed in the Appendix. There 
are also solutions in which the orbital distribution in the 
mean field configuration of the ground state is not unique. 
For instance, among the degenerate states at large 7 we 
find the polar vortex with (nio) — 1 an d a configuration 
with (n 20 ) = (n Q0 ) = \. 
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E. Comparison with finite size exact states 



It is of interest to try to match the mean field states 
with ground states found in diagonalization studies, such 
as those shown in Figure[21for the disc. Since these have 
definite values of the quantum numbers, they can be com- 
pared with the mean field states only by projecting the 
latter to components with definite quantum numbers 33 . 
For each N and L, the value of the spin picked out should 
reflect the form of the interaction, and should presumably 
be the maximal value in the ferromagnetic regime, and 
the minimal value in the antiferromagnetic. For the spin- 
independent case c 2 = 0, the lowest SU(3) ap i n quantum 
numbers are favored as ground states. 

We will not attempt to identify all the states in Fig. [21m 
this way, but only some of the more prominent. We have 
already mentioned that the mean field state at I = 1 and 
c 2 — corresponds to the BTC singlet state. Since the 
mean field state has equal mean occupation of (to, a) = 
(0,|), (1,0), and (2, f), it does contain a unique singlet 
component which is exactly the BTC state. When c 2 
is turned on, the quantum numbers remain at (L, S) — 
(6,0), but the state will be slightly altered in its details. 
The corresponding skyrmion spin textures on the sphere 
are also discussed in the Appendix. 

The BDC multiplet at L = N/2 for c 2 = that uses 
only to = 0, 1 also deserves comment. This corresponds 
in mean field theory to the N v = 1 case discussed above 
and in the Appendix. When c 2 < 0, it becomes a half- 
skyrmion or meron, which survives for all — 1 < 7 < 0. 
This meron has no projection to spin 0, and anyway for 
this regime maximal spin is expected in the ground state. 
Indeed, for N — 6 the corresponding (L — 3) state has 
S = 3. The whole regime L/N < 1 for ferromagnetic 
interactions resembles what one expects for skyrmions, 
that is L (corresponding to N — L on the sphere) de- 
creasing as S increases, as S = N — L^. For c 2 > and 
L = N/2, the lowest-spin part of the BDC state becomes 
the ground state. 

At larger positive 7, there is a prominent region of 
(L,S) — (6,0) in the N — 6 data. At the largest 7, we 
expect that this can be identified (in the same sense as 
the preceding discussion, or as in RefJ&) with the polar 
vortex state of this section. (In a finite size study, one 
would not expect to see a transition from the BTC state 
at c 2 — to this polar vortex with the same quantum 
numbers at large 7.) The jump from L = to L = 6 
expected from the mean field is seen in Fig. [5] At smaller 
7, a (3, 0) region is seen. We speculate that this state 
corresponds to a single 7r-disclination at the center of the 
trap (with a second one at infinity, or the opposite pole on 
the sphere), and that the region corresponds to i\ < £ < 
2/3 in the mean field results. Notice also the prominent 
semi-plateaus near I = 0.5 in the plots in Fig. at larger 
7- 



VI. VORTEX AND SKYRMION LATTICES 

Upon driving the system faster, multiple skyrmions are 
induced. These are expected to form a lattice and can 
be well treated in a (quantum Hall) mean field analysis. 
Such an analysis was performed by Kita et alH, who 
found a range of different lattices for c 2 > 0, depend- 
ing on the relative strength to cq and the rotation. By 
including higher Landau levels, they were able to show 
that some of these lattices are qualitatively identical at 
high and low rotation. Near c 2 = 0, the (scalar) vortex 
breaks up into three vortices, one for each spin compo- 
nent, forming a triangular lattice. For c 2 > 0.069co, the 
vortex splits into two 7r-disclinations, which make up a 
square (anti-ferromagnetic) lattice. 

We have carried out a program, similar to Mueller and 
Ho2&, appropriate for a mean field LLL description of a 
multi-component condensate. The LLL approximation 
(in the limit u> — > u>q) fixes the vortex lattice spacing to 
be equal to the harmonic oscillator length. Note that 
this is different from the Thomas-Fermi regime, where 
the distance is fixed by the number of vortices, as the 
density of bosons is the same as in a non-rotating trap. 

Under the assumption that the vortices in each spin 
component form a Bravais lattice, we can choose the one- 
particle wavefunctions to be the torus wavefunctions with 
N v flux quanta (typically, N v = 1 or N v = 2). For a 
scalar condensate, the lattice is completely specified by 
the geometry r of the torus. This wavefunction is peri- 
odic up to a gauge transformation, equivalent to requiring 

A'(r') =A(r). 

In the case of multi-component condensates, however, 
more general boundary conditions are possible. We only 
need to demand 

^'(r + Li) - e iA *W[/^(r), (21) 

where Lj (i = 1, 2) define the geometry, and Aj is the 
gauge transformation mentioned above. The matrices U\ 
and £/ 2 should commute, U\UiU\ = 1) to obtain 

single- valued wavefunctions. 

We require that U± and t/ 2 commute with the Hamilto- 
nian, so that the energy of a unit cell is well-defined. For 
7^0, this implies f/, € SO(3) sp i n . The common eigen- 
vectors of U\, U2 then have eigenvalues (1, e llpi , e~ t(pi ) 
and (l,e^ 3 ,e-^ 2 ). With an overall SO (3) rotation, we 
can fix the direction of the vector with eigenvalue 1 to be 
parallel to z in spin space. With this, the unit cell of the 
magnetic order (seen in the spin density which is gauge 
invariant, for example) is larger than that of the density, 
but always contains an integer number of the latter. 

Using this approach, we can confirm a large part of the 
phase diagram of Kita et a&&, but we also find additional 
phases in the ground states at large 7. These are polar 
phases, for which we use a unit cell with a single flux 
quantum. We will use <f> — arctan7 as the parameter. 
The minimization procedure uses a simplex downhill al- 
gorithm in the geometry r and the phases ip±, </? 2 . The 
wavefunction is obtained from the polynomial free energy 
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FIG. 10: The different lattices found in rotating spin-1 boson condensates. The first picture is the mapping of the spin-sphere 
to the colors used. Top and bottom correspond respectively to the north and south pole. The intensity codes 1(5)1, the size of 
the spin-vector. Other pictures are the spin expectations at different ratios cijcc,: <j> — —0.1, —0.05, 0.01, 0.016, 0.04, 0.1, 0.54, 
0.7 and 0.9. The last picture shows the density, as the spin vanishes. 



by using a conjugate gradient algorithm, starting from a 
random point. The wavefunction in general is unique up 
to a phase and a SO(3) rotation along the z-axis. 

The phases we obtain, as illustrated in figure El are 
as follows: 

ferro lattice. A major part of the ferromagnetic phase 
diagram (— 7r/4 < tj> < —0.08) is covered by a lat- 
tice with N v = 2 flux quanta in the unit cell. This 
is the same lattice as one obtains for the spin-^ 
bosons with full SU(2) symmetry or, equivalcntly, 
the quantum Hall ferromagnet with the Lande fac- 
tor <7 = 0. If we consider the spin-1 to be com- 
posed of two spin- 1/2 particles, then N v = 2 for 
the spin-1 bosons corresponds to N v = 1 for the 
spin-i particles. This structure is related to the 
N v = 2 skyrmions discussed in the Appendix. 

skyrmion-vortex lattice. At <fi —0.08, it becomes 
beneficial to include vortices ("merons"). The unit 
cell now has N v = 3, with both a skyrmion and a 
vortex. Based on direct computations in disc geom- 
etry (see below) , we expect that this phase does not 
extend to cf> = 0, but that there are other phases in 
the weakly ferro regime —0.02 < <f> < 0. 

triangular vortex lattice. Exactly at cf> — 0, the nodes 
in the three components are arranged in a triangu- 
lar lattice. This lattice can be realized with N v = 3 
and (fi — (f2 — 0. The mean-field components b mfll 
(m = 0,1,2), form a unitary (C/(3)) matrix. This 
lattice is not shown in figure ITU1 as the S'0(3)-spin 
is not well defined. 

square ladder. The triangular vortex lattice of the C2 = 
case is essentially unchanged up to cj) = 0.0143, 



being squeezed only. However, the SU(3) symme- 
try is broken. This spin shows a ladder structure, 
where adjacent ladders are shifted by 3/2 rung- 
spacings. 

canted ladder, 0.0143 < (j) < 0.0193. The ladder struc- 
ture stays intact, however the rungs are now canted. 

triangular ladder, 0.0193 < < 0.069. 

square 7r-disclination, at <f) w 0.069, there is a first 
order phase transition to the square 7r-disclination 
lattice. Only the | and J, components are present 
in this lattice. 

squeezed 7r-discl., 0.428 < < 0.62. The lattice 
is squeezed in one direction and expanded in the 
other. 

triangular 7r-discl., 0.62 < 4> < 0.786. At w 0.62, 

there is a first order phase transition to a triangular 
7r-disclination lattice. 

polar Abrikosov, beyond cf> « 0.786, the 7r-dis- 
clinations are unstable and the systems prefers to 
have only one component, such that (S) = ev- 
erywhere. The vortices of this component form 
an Abrikosov lattice, with vanishing density at the 



The phases at (j> > 0.428, and at </> < have not been 
observed before. Figure ITU1 shows the spin texture in the 
various lattices, with colors coding the direction of the 
spin vector and the intensity marking its length, so that 
black regions indicate places where all components of the 
spin vector vanish. [For the lattice at <fi = 0.9, which 
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is the polar Abrikosov lattice, the spin density vanishes 
and we plotted the particle density instead.] The particle 
density is finite in all lattices except the polar Abrikosov 
one. 

To check whether the Ansatz is sufficiently general in 
the complete phase diagram, we have supplemented the 
above analysis by direct numerical computations of LLL 
mean field ground states in a disc geometry, with lu < 
loq. Since no periodic structure is imposed, the lattices 
form spontaneously. These computations show that the 
torus correctly reproduces the dominant phases such as 
the square lattice of 7r-disclinations and the skyrmion and 
skyrmion-vortex lattices. In the region —0.02 < <f> < 
the two geometries showed different lattice structures, 
possibly due to finite size effects. We leave conclusive 
results in this region for future work. 

At special values of ipi, (p 2 , when they are both of 
the form pir/q (p, q integer), it is possible to realize the 
lattice by using a larger unit cell and identical phases 
for all three spin components. An example of this is the 
triangular lattice at c 2 = 0, where <p\ = —(p 2 = 2tt/3. 
In this case, we can realize the same lattice by using a 
torus with 3 flux quanta and ipi = ip 2 — 0. The other 
example is the square 7r-disclination lattice, which can be 
described by using 2 flux quanta. This can be compared 
to the spin-1/2 situation^, where the lattice at gi 2 = 
.9i = .92 (unbroken SU(2) symmetry) can equivalently be 
described by using a torus with 2 flux quantai2i2£. 



VII. QUANTUM HALL LIQUIDS 

In the LLL approximation, the filling fraction v de- 
fines the average number of bosons in an orbital. Upon 
increasing the rotation further, and thus reducing the 
discreteness of the occupation numbers becomes impor- 
tant. Mean field theory becomes less useful due to quan- 
tum fluctuations. At some point, the condensate is de- 
stroyed, and a sequence of quantum fluids takes over. In 
the scalar case, the critical v was estimated to be v c w 10 
by the Lindemann criterion^ (that is, the average fluc- 
tuation in the position of the vortices equals the separa- 
tion between them.) Explicit calculation of small systems 
have confirmed this transition and found v c w 6— 10. A 
similar transition will occur in the spin-1 case, although 
we have not calculated the appropriate v c . We expect 
this to be of the same order of magnitude as in the scalar 
case. In the present section we investigate the quantum 
Hall fluids that appear within region II of the phase dia- 
gram in this regime. 

In the extreme limit u> — > ujq (is — » 0), we can analyze 
the quantum liquids analytically in a part of the phase 
diagram, as we can explicitly find the zero-energy eigen- 
states of the Hamiltonian. Two of the series we propose, 
the SU(4)k and the SO(5)k series, have a member of this 
form for k = 1. The third series consists of a generaliza- 
tion of a family of fractional quantum Hall (QH) states, 
the hierarchy /composite fermion states, to spin-1 parti- 



cles. We present some numerical results on small sizes, 
which unfortunately are probably not conclusive for the 
nature of the states, due to the restriction to insufficiently 
large sizes. 

A. SU(4)k series 

It is straightforward to construct zero-energy states at 
c 2 = 0. The repulsive contact interaction dictates that 
the wavefunction should have a node whenever 2 particles 
are on the same place. Furthermore, two bosons with the 
same spin should have a double zero in order to maintain 
a symmetric wavefunction. In terms of the components of 
the wavefunction for given values fi = x, y, z of the spin 
for each particle, we can write such a state down^L^ (as 
in Refpiii, we add a tilde to the wavefunction to indicate 
that it has to be multiplied by the usual Gaussian factors 
for the plane, or the rational factors for the sphere — see 
Section ITTft : 

^2,2,2,1,1,1/ x x v V . z 2 \ _. 

n n#-*?> a n ik -*?"). ^ 

where n a denotes the number of particles with spin a. 
The lowest angular momentum L for which this state can 
be realized is when N x = N y = N z = N/3. Notice that 
L = N(2N/3 — 1), and the filling factor, which can be 
defined as v = limjv^oo N 2 /(2L), is v = 3/4. This state 
is a straightforward generalization of both the Laugh- 
lin v — 1/2 state for spinless particles, and the (2,2, 1) 
Halperin spin-singlet state for spin-1/2 bosons. The full 
state for N x — N y = N z — N/3 is an SU(3) sp i n sin- 
glet ((p, q) = (0, 0)), and arguments similar to Laughlin's 
plasma mapping show that moreover it has short range 
spin correlations. Because the state vanishes whenever 
the particles are at the same point, this state is a zero- 
energy eigenstate for all values of Co, ci- However, it will 
only be the ground state (at L = N(2N/3 — 1)) when 
-ffint is positive, that is in the region go, g 2 > 0, within 
region II. For larger values of L, there are many more 
zero-energy eigenstates, so the ground states are degen- 
erate in the window go, g 2 > 0. On minimizing H u with 
respect to L, this implies that the lowest possible filling 
factor as u> — * luq from below is v = 3/4 in this regime 
(within the model in Sec. ILT|) . 

The state in eq. (|22|l is an exact eigenstate, but in gen- 
eral we are not able to find the exact highly-correlated 
ground states of H lnt . Instead, we seek to understand nu- 
merical results, and make predictions for the physics at 
larger sizes by using (among other techniques) trial wave- 
functions. These states, which are not generally exact 
for any know two-body interaction, serve as paradigms 
for the phases of matter in the thermodynamic limit, as 
they possess interesting ("universal") properties such as 
the quantum numbers and statistics of their excitations 
that are robust against small changes in the Hamiltonian, 
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until some phase boundary is passed (this philosophy has 
been discussed e.g. in Refii^). One way to produce trial 
wavef unctions, that is closely connected to their univer- 
sal properties, is to use conformal field theory (CFT). 
We will show how to obtain wavefunctions from a CFT 
in somewhat more detail in the next section. The CFT 
describing eq. 1221) is SU(A)i, so following a strategy in 
Read-Rezayii^, and motivated by the analogous results 
for scalar bosons^, we can consider a series, SU{A)k, 
where k — 1, 2, .... These states have filling factor 
> 3/4, and are explicitly SU(3) sp - m invariant. Hence we 
expect them to be relevant near ci = 0. 

The trial states we consider have wavefunctions, in spin 
components (generalizing those in22iiii to spin-1), 



*({*}) = S ( 



roups 



n 



2,2,2,1,1,1 



(23) 



In this construction the N = 3pk bosons are first parti- 
tioned into k groups, each with p particles of each spin 
component x, y, z. 

For each group we write a Halperin vp 2 ! 2 ' 2 ' 1 ' 1 * 1 fac- 
tor, and these are multiplied together. Finally, the sym- 
metrization operation S over all ways of dividing the par- 
ticles into groups is applied. The angular momentum is 
L = N[2N/(3k) - 1], and the filling factor (as N -> oo 
at fixed k) is therefore v = 3k/ 4. It happens that if we 
put k = N/3, the state contains k groups of three bosons 
each, and the state is exactly the L = N BTC. How- 
ever, we do not believe this is particularly significant, as 
the BTC state is the unique SU(3) sp i n singlet state at 
L = N. 

The states (|23|l are zero-energy eigenstates of a Hamil- 
tonian consisting of a k + 1-body interaction: 



H su(i) k - V 



il<---<ik+l 



(24) 

The interaction is positive for V > 0, so all eigenstates 
have E > 0. This interaction penalizes k + 1 particles 
at the same point. Therefore, zero-energy eigenstates 
are those in which the wavefunction vanishes if any k + 1 
coordinates coincide, regardless of the spins. One can see 
that the above function has this property, even before the 
symmetrization over groups, as for any k + 1 particles, at 
least two must be in the same group, forcing the function 
to vanish. For less than k+1 particles at the same point, 
it does not necessarily vanish. In fact, for each k, (|23|) is 
the unique zero-energy eigenstate of Hgu^ k with lowest 
angular momentum L. 

For the same Hamiltonian on a torus, there are again 
zero energy states, at least for N divisible by 3k. For 
these cases, the degeneracy of these SU (4)^ ground states 
is 



# fc = i(fc + l)(fc + 2)(fc + 3) 
6 



(25) 



We have verified that this result, which can be inferred 
from the CFT connection, is reproduced by exact diago- 
nalization of the Hamiltonian eq. I|24|) on the torus. 



Like other incompressible QH states, the phases of 
matter exemplified by the trial states l|23|l possess point- 
like quasiparticle excitations which may have fractional 
particle number (relative to the background density) 
and/or spin. The particle number associated with the 
elementary quasiparticles can be found once it is under- 
stood that, similar to the RR states 1 ^, the SU(A)k states 
are clustered states, in which particles occur in clusters 
of 3k (in an 5t/(3) sp i n -singlet). Then a similar argu- 
ment to that given in Sec.[I]shows that they carry charge 
q = ±1/4. They also have spin 1. The quasiholes, at 
which there is a deficiency of particle number, can be 
studied as zero-energy eigenstates of Hguu\ h , and fairly 
explicit trial wavefunctions can be found using the rela- 
tion with CFT. 

The statistics of quasiparticles in 2D can be defined 
in terms of adiabatically dragging them along paths, 
keeping them well-separated, to exchange them. For 
the SU(4)k states, in the case fe = 1, the statistics 
are "Abelian" ; the wavefunction acquires a phase factor 
when two particles are exchanged, just as for the Laugh- 
lin states. For fc > 1, however, the statistics becomes 
"non-Abelian"^*ii. This means that there is a degener- 
acy, when the positions of the quasiparticles are fixed (in 
general, this is true only when the quasiparticles are well- 
separated, though for the k + 1-body Hamiltonian and for 
quasiholes it is exact for any separation). In terms of the 
trial functions (|23|) and their generalization to quasiholes, 
this degeneracy is caused by the symmetrization proce- 
dure S, which destroys the "group" quantum number of 
the (quasi-) particles 39 . The general framework to ob- 
tain these degeneracies from CFT has been worked out 
and agrees with numerical results for the Moore-Read, 
RR, and other stateo 17 i 41 i 42 i 43 i 44 . We expect the same 
framework to apply here. When the quasiparticles are 
exchanged adiabatically, the effect is a matrix operation 
within these spaces of degenerate states, described by the 
braiding matrices of the corresponding CFT— — hence 
the term non- Abelian statistics. 



B. SO(5) k series 

Similar to the SU(4)k case, the SO(5)k states can be 
written in the form: 



50(5) 



where now 



^ffi (5) ({^}) = Pf 



({*}) = $ 



groups 



-(groups w fc = i 



(26) 



• o 



Hizi-Zj). (27) 



i<j 



Here the spin states for all the particles are included ex- 
plicitly (the product of spin states Q being the tensor 
product) Pf Mij denotes the Pfaffian of an antisymmet- 
ric matrix . In the present case the N = 2kp particles 
are partitioned into k groups, with 2p particles in each. 



B 50(5) fe series 

The particles in each group form an 50(3) sp i n singlet. 
The product over these groups is then symmetrized. The 
state as a whole is clearly an 50(3) sp i n singlet, and has 
angular momentum L = N[N/(2k) — 1], so the filling 
factor is v = k. 

The k — 1 case, which closely resembles the Moore- 
Read paired stateAi but for spin-1 particles, is the ex- 
act ground state of our two-body Hamiltonian H lrA when 
<?o = 0. That is because, in the state vj)^ ^^ two parti- 
cles are found at the same point only if they have total 
spin 0. Again, the ground state as above is the unique 
zero-energy eigenstate at the stated angular momentum, 
but at larger L there are many more zero-energy states. 
So v — 1 is the lowest filling factor possible at go = 0. 
This implies that in finite size on the disc a boundary 
between ground states with the L values of the 5t/(4)i 
and 5*0(5)1 states must run into ui = uio at g = (this 
is for N divisible by 6, but there will be similar state- 
ments for other values). Such behavior is seen in Fig. El 
For go < (7 > 1/2), we do not know the exact lowest v 
that occurs as uo — > loq from below. The 50 (5) 1 state can 
be interpreted in terms of BCS spin-singlet complex-p- 
wave pairing of composite fermions, in which the Pfaffian 
represents the pairing in position spac o 17 i 20 . 

More generally, for each k there is a Hamiltonian for 
which the 50 (5) & states are exact zero-energy eigen- 
states, again given by a k + 1-body interaction: 

H so(5) k = V S(z n - z i2 ) ■ ■ ■ S(z ik -z ik+x ) 

il <•••<«*. +1 

xP fe+ i(ii,...,i fe+1 ) . (28) 

This interaction includes a projector Pk+i(ii, . . . ,ik+i) 
of the spin state of the k + 1 particles concerned onto 
total spin k + 1. 

For general k, these states can be considered to be 
built up out of clusters of 2k particles in a spin singlet. 
From this fact we can obtain the fractional particle num- 
ber of the elementary quasiparticles, q = ±1/2. Also, 
the quasiparticle spin is 1/2, which is fractionalizcd com- 
pared with the spin 1 of the underlying bosons, and so 
the number of these quasiparticles must be even. For 
k = 1, there are also excitations with zero particle num- 
ber that behave as fermions with spin 1. In this case, 
the universal properties may be understood by a simple 
extension of the methods of Refill to this case. 

Because it is difficult to see through the symmetriza- 
tion operation 5, we will provide some details on the 
conformal field theory behind these states. Such a CFT 
description allows us to obtain more insight into the topo- 
logical properties, such as degeneracies and braiding. For 
example, to obtain the degeneracy of ground states on 
the torus, CFT tells us that we only need to know the 
number of non-trivial representations at level k. In the 
case at hand, this number turns out to be |(fc + l)(fc + 2). 
Again, we verified this number using exact diagonaliza- 
tion of the (k + l)-body interaction. 

The chiral algebra of the CFT which describes these 
states is based on the 50(5)& affine Kac- Moody algebra. 
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FIG. 11: Roots and weights of the algebra 5*0(5). The con- 
densate operators are associated to roots (filled symbols) 
and the fundamental excitations (f> correspond to the weights 
of the spinor representation (open symbols) . 

50(5) is a rank 2 Lie algebra, which contains mutually 
commuting 50(3) and U(\) Lie subalgebras, which we 
can identify with the symmetries under 50(3) sp i n and 
number conservation. In these subalgebras, the genera- 
tors are respectively 5 + , 5~, S z and c respectively. Ac- 
cording to the CFT-qH correspondence^, we can also 
obtain the quantum Hall state wavefunctions as corre- 
lators in the chiral part of a CFT, in which the parti- 
cles (bosons) should be represented by fields that have 
Abelian braiding properties. In the present case (and 
the SU (4)fc case is similar), the bosons correspond simply 
to a different triplet of current operators of the 50(5)fe 
affine Kac-Moody algebra. Thus the wavefunction, now 
in spin components, can be written 

*({*}) = lim z^ k (e- lN ^^( Zoo ) 

Zoo— » OO 

J ai (zi) ■ ■ ■ J aN (z N )) 1 (29) 

with J a (a —I, 0, |) an 50(3) sp i n triplet of currents in 
the affine Lie algebra, that carry t/(l) charge +1. The 
currents are shown in the 50(5) weight lattice (fig. fTTfl . 
as ^i- The operator whose position tends to 

00 represents a background charge, such that the total 
U(l) charge of the operators is zero (as it must be in 
order that the correlator be nonzero). The currents can 
be expressed as 

J a (z) = ^(zK^/^M, (30) 

where ip a is a parafermion field, of conformal weight 

1 — 1/fc for the long roots and 1 — 1/2/c for the short 
roots. The vertex operator e l/3 ° "$l^(z) contains the two 
free boson fields ip c (charge) and ip s (spin), with (3 the 
position in the root lattice. For J a these are /3| = (1, 1), 

A) = (1,0), 0i = (1,-1). 

The parafermions simplify when we specialize to the 
case k = 1, where they reduce to the identity operator for 
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the long roots (vPf and Wj, in particular) and a Majorana 
fcrmion for the short roots (^o)- The correlator can then 
be readily written down, as the correlation functions for 
Majorana fermions is well known: 

(#i)-^W) = Pf (—?—)■ (31) 

V z i *j / 

This reproduces the SO(5)i wavefunction, in the same 
way as for the spinless Moore-Read stated. 

We note that the same 50(5)^ algebra was used in 
another construction^, which was for a system of spin- 
1 /2 particles. The present case differs in that the physical 
50(3) sp i n symmetry is embedded differently in 50(5), 
because of the different spin of the underlying particles. 

Wavefunctions for zero-energy states containing quasi- 
holes can also be written down as chiral correlators, 
which now contain vertex operators for primary fields 
of the SO(5)k algebra that represent the quasiholes. For 
the k — 1 case, these contain (in the scalar field plus 
Majorana fcrmion language) a spin field for the Majo- 
rana fermion, and give rise to quasihole wavefunctions 
analogous to those for the Moore-Read stat o 1 ?' 41 . 



C. Composite fermions 

Alternative QH states to the rather exotic series in the 
previous two subsections can be constructed by apply- 
ing conventional methods to spin-1 bosons. One such 
approach, as in the case of scalar bosons, is to map the 
bosons onto composite fermions, by attaching (say) one 
vortex to each boson. These fermions see a reduced ef- 
fective magnetic field, and one can construct an incom- 
pressible state when an integer number of Landau lev- 
els in the effective magnetic field is filled with all three 
components. This construction gives states with filling 
factors v = 3p/(3p± 1), which are SU(3) spm singlets. 

One can interpret the Moore-Read state as the p-wave 
pairing of composite fermionsit In this case, p-wave 
50(3)-singlet pairing is possible (in contrast to the spin- 
1/2 case) and indeed, we have seen that the 5*0(5)1 state 
can be interpreted this way. In the SU (3) symmetric case 
at v = 1, no 2-particle 5C/(3)-singlet pairing is possible 
and there are two options for the system. One is to form 
a Fermi liquid, the other is to spontaneously break the 
symmetry and form [p, q) = (2, 0) pairs. Note that this 
last possibility includes the 5*0(3) singlet and thus can 
be continuously connected to the 50(5)i state. 



D. Vortex lattices without polar order, and 
nematic QH liquids 

The earlier discussion of QH liquid states focused on 
singlets under SO (3) sp i n , with short range spin correla- 
tions. It is interesting to wonder also if QH liquids with 
some form of spin ordering could occur. One possibility 



would be a ferromagnetic QH liquid. Such states can be 
easily written down, by using any wavefunction for a QH 
state of spinless bosons, with all the boson spins in the 
a =t state (or a global spin rotation of this). One might 
expect these to occur in the ferromagnetic (02 < 0) part 
of region II, but in fact we see no sign of them: leav- 
ing aside the skyrmion textures in the BEC at low L, 
at larger L all ground states are spin singlets. We note 
that for spin-1/2 electrons, spin-polarized states can oc- 
cur, e.g. at v = 1, even for spin-independent interaction, 
due to exchange effects. However, the exchange effects 
are presumably different for bosons. 

A more feasible-looking possibility is QH states with 
polar spin order, perhaps in the antiferromagnetic region 
C2 > 0. In the regime at large 7 where mean-field the- 
ory predicts the Abrikosov vortex lattice, the spin-order 
is polar. In the polar state, the vector condensate can 
be written as (Vv) = e%lpn iJ.i with tp the phase and n a 
real vector. In Abrikosov lattice, the magnitude of the 
vector h and the phase ip vary to give a triangular lattice 
of vortices. We can now imagine that quantum fluctu- 
ations destroy either part of the order (restoring cither 
the phase or the spin-rotation symmetry) without the 
other. When the U(l) and translational symmetry that 
are broken in the vortex lattice are restored, the ground 
state is a QH fluid. For large quantum fluctuations one 
might expect that the QH liquid has restored 5J7(2) spin 
symmetry. However, the two transitions at which these 
symmetries are restored are independent, and the tran- 
sitions could in principle occur in either sequence as we 
go to smaller v. The intermediate phase in which spin 
symmetry is restored but not the phase would be a vor- 
tex lattice in a boson paired state, a BEC of boson pairs. 
This would be characterized by having a nonzero expec- 
tation value of VvWVvW' which is invariant under 
%jj — > —ip. (Such a vortex lattice would also be possible 
with vortices containing a half-unit of vorticity each, in- 
stead of integers as we assume otherwise, and this might 
be reached by restoring symmetry in the 7r-disclination 
lattice state.) 

The other possible sequence of transitions would be 
that in which the intermediate phase is a QH liquid with 
restored translational and phase symmetry, but still has 
the polar order, which breaks SU(2) spin . The single- 
boson expectation (ipfi(z)) would be zero, but if we look 
at the composite operator ipj 1 (z)ip^(z), this can have an 
expectation value. This matrix has a trace equal to the 
density, which is uniform by assumption. The trace- 
less Hermitian matrix obtained by subtracting off the 
trace contains antisymmetric and symmetric parts. The 
(imaginary) antisymmetric part corresponds to a spin-1 
irreducible tensor that is simply the spin density, which 
is assumed to be zero here. The (real) symmetric part 
corresponds to a spin-2 irreducible tensor. This is the 
order parameter of a polar or nematic state, which rep- 
resents a vector h, but is invariant under n — > —h, so 
it parametrizes 5 2 /Z2 = RP 2 . Trial wavefunctions for 
these nematic quantum Hall states can be written down 
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as those for scalar bosons, times a spin state such as 
a = for all bosons, or as a spin-rotation of this. 

It would not be surprising if such nematic QH liq- 
uids occurred in the phase diagram at large 7, now that 
the corresponding (polar Abrikosov) vortex lattices are 
known to be present. In finite size, the ground state 
would always be low spinS, S = or 1, and there need 
be no transition separating it from a state at the same L, 
S with short-range correlations, such as the SO(5)i state 
at go = 0. Thus the appearance of such nematic order in 
a QH fluid of spin-1 bosons in the thermodynamic limit 
cannot yet be ruled out. 



E. Numerical results 

To examine how well the proposed states describe the 
true ground states, we have performed exact diagonal- 
ization of small systems. In the regime L < N, we have 
used both the disc and sphere geometries. As we have 
seen, these results differ somewhat. But when looking at 
fast rotation, however, where the filling factor is of or- 
der 1, the system is spread out into a pancake. It makes 
sense to focus attention on the interior of the disc and 
avoid edge effects. This can be done by using an edgeless 
geometry such as the sphere or torus. Here we will be 
interested in the ground states in which (unlike the work 
earlier in this paper) we find the ground states without 
constraining L. QH liquid ground states will usually then 
show up as L = states. At finite sizes on the sphere, 
such ground states that form a sequence of sizes tending 
to a particular filling factor v in the thermodynamic limit 
lie on a sequence of the form N v = N/v — s2£ . Here s 
is known as the shift, and its appearance is connected 
with the coupling of the particles to the curvature of the 
sphere. The value of s depends on the liquid state, not 
only on v. N v can be obtained from the angular momen- 
tum on the disc as L = NN v /2 (all states have L = 0). 
That is, when the states in the \& notation are written for 
the sphere, one can take N v as small as possible, so that 
N v equals the highest power of any Zi appearing in the 
wavefunction (see Sec. [HJ|. For QH ground states, this 
will usually mean that L = 0. For the SU{A)k ground 
states, we have N v = 4N/(3k) - 2, while for the SO(5) k 
ground states N v = N/k — 2. We will compare the results 
of numerical solution for the ground states with these se- 
ries of trial states. 

The SU(A)k states with k = N/3 are the exact ground 
states for c 2 — at N v = 2 on the sphere. It is not 
surprising that they are eigenstates, because they are the 
only SU(3) singlet states, as in the case of L = N for the 
disc. 

As a further test for the SU(4:)k states, we have looked 
at sizes (N,N V ) at which such a ground state could lie 
for k > 1. Since N must be divisible by 3k, such sizes 
increase rapidly even for k = 2. The next case after the 
trivial BTC for k = 2 is N = 12, N v = 6. Here it turns 
out that the overlap-squared of the exact ground state 



for c 2 = with the trial state is&. \(SU(4) 2 \G.S.)\ 2 = 
0.915226. 

The SO (5) 1 state was shown to be an exact ground 
state for cq > 0, go = 0. The higher members of this 
series, however, had a vanishing overlap with the ground 
states throughout the phase diagram. 

As a further test of our proposed wavefunctions, we 
have performed calculations for torus geometries. On 
the torus, one simply has v = N/N v for the finite size 
sequence of ground states that tend to a fluid of filling 
factor v in the thermodynamic limit. We saw in the mean 
field analysis of the skyrmion lattice, that the lattice can 
only be observed when the number of flux quanta is a 
multiple of three. However, at low filling factors, we ex- 
pect to see quantum Hall states at v = 3fc/4. To be 
able to observe these, N v has to be a multiple of 4. Un- 
fortunately, this implies torus sizes which are too large 
to observe both the quantum liquids and the skyrmion 
lattice. 

To see if the proposed wave-functions are good can- 
didates, we are therefore forced to look at tori which 
frustrate the mean field skyrmion lattice. The cases we 
considered are N v = 3, 4, 6. For N v = 4, we find that 
the ground states are exactly given by the SU{4)k series. 
However, as for the BTC states on the sphere, this is 
due to the fact that the trial ground states, which are 
5J7(3) sp in singlet states, span the space of all SU{3) sv - m 
singlets on the torus, which has dimension equal to 
eq. (|25|) . the degeneracy of SU(A)k torus ground states. 
Clearly this must be independent of the geometry of the 
torus (described by r), and we verified this in some cases. 

In fig. ^] we have plotted the particle-hole excitation 
gap 

A { N)=N( E ^ N - 1 K E{N+1) -2^1), (32) 
1 ' V N - 1 N + l N J 

where E(N) is the ground state energy for N particles, 
for N v = 4. In the thermodynamic limit, this quantity 
will exhibit upward peaks at filling factors that corre- 
spond to incompressible states. We performed the diag- 
onalization for total SO(3) spin S = 0, 1, 2. 

For N v — 6 ( figure IT3)) . we focussed on the state at 
v = 3/2, N = 9, which corresponds to k = 2. We have 
calculated the overlap-squared with the SU{4) 2 state to 
be 0.939804. Another feature in the N v = 6 plot is the 
state at v — 1 (N — 6). This could be a precursor to 
a paired composite fermion state; however, the overlap 
with the SO state was small. 

F. The boundaries of region II 

The behavior at the phase boundaries I-A/II and I- 
B/II (see figure ^1 deserves special attention. At the 
boundary I-B/II, where 7 — c^jc^ = —1, the Hamilto- 
nian simplifies as g 2 — and only the contact interaction 
which projects onto the spin-singlet channel remains. As 
a result of this, large degeneracies occur. For example, 
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FIG. 12: Particle-hole excitation gap A(iV) versus u, for 
iV„ = 4, in a rectangular geometry, a/b = V3/2. The peaks 
can be interpreted as an indication of incompressibility of the 
corresponding states. For N = 3,6,9 we verified that the 
ground states, which are degenerate, are exactly the SU(4i)k 
quantum Hall trial states with k = 1, 2, 3, respectively. 




FIG. 13: Particle-hole excitation gap A(JV) versus u, for N v = 
6, in a rectangular geometry, a/b = v3/3. 



all fully-polarized states (S — N) have zero energy. We 
have not obtained analytic expressions for these degen- 
eracies. That they are not due to the specific geometry 
was observed on the torus. The zero energy states are 
not sensitive to changes in the geometry. As examples of 
these degeneracies, we have in figure ITU tabulated the L, 
S quantum numbers of the zero energy states for N = 6, 
N v = 2 and for N — 5, N v = 3 on the sphere, and for 
N = 6 particles in the disc geometry. 
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FIG. 14: Degeneracies of zero-energy ground states on the 
sphere at 7 = -1 for N = 6, N v = 2 (top left), JV = 5, 
N v — 3 (top right) and N = 6, N v = 00. All multiplicities 
refer to highest weight states of the orbital 50(3) symmetry. 



(i) low rotation, such that the angular momentum L is 
less than or equal to the particle number N, where the 
system is beginning to contain some vorticity; (ii) higher 
rotation, where the bulk of the fluid accommodates vor- 
ticity and is occupied by a lattice of (possibly coreless) 
vortices, which we considered as infinite periodic struc- 
tures; (iii) the quantum Hall regime, in which the vortex 
lattices are replaced by translationally-invariant quan- 
tum fluids, which we considered in edgeless geometries 
(unfortunately, for spin-1 bosons the finite size restric- 
tions are here very severe). The transition regions be- 
tween these regimes, namely that in which the system 
contains a small number (larger than two) of vortices, 
and that at the critical filling factor at which the vor- 
tex lattices are replaced by the quantum liquids, were 
not considered. The results show a rich variety of phases 
as the interaction parameters, especially the ratio of the 
coefficients of spin-dependent and spin-independent in- 
teraction terms, are varied. The results obtained here, 
especially those at lower rotations which should be more 
easily accessible, should motivate further experiments to 
rotate spin-1 bosons with unbroken spin-rotation sym- 
metry. 



VIII. CONCLUSION 

In this paper we have studied the phase diagram of 
spin-1 bosons in a rotating trap, within the LLL approx- 
imation, using a variety of techniques (numerical diag- 
onalization, mean field theory, and analytical construc- 
tions). We concentrated on certain regimes. These were 
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APPENDIX A: CLASSIFICATION OF 
TOPOLOGICAL DEFECTS 

Here we continue the discussion of the topological clas- 
sification of defects or excitations. 

In Sec. [i] we explained the appearance of two types of 
ordered BECs for spin-1 bosons. In these the order is 
constant in space. More generally, a condensate will pre- 
fer to have the same type of order locally at almost all 
points in order to lower the energy, but there are types 
of excitations in which the order as described by a point 
in the order parameter manifold can vary in space. Ex- 
citations of this type that are stable under continuous 
deformations of the order (known as topological defects) 
can be classified by methods from topology. One type 
is those in which the order breaks down (possible, the 
density goes to zero) at a point in space (we consider two 
space dimensions here) . These are classified by the funda- 
mental or first homotopy group 7Ti of the manifold. Such 
objects exist for example in the case of a one-component 
(scalar) condensate, where the order parameter manifold 
is the circle S 1 . For vector condensates, an example is 
the polar state, for which 7Ti(5 1 x S 2 /Z2) = Z, the group 
of integers. Because these involve the phase of the con- 
densate (related to the S 1 part of the manifold) winding 
as one encircles the singular point, these objects in both 
examples are vortices that are relevant when vorticity is 
forced into the system. In the polar state, the vortices of 
the smallest vorticity carry a half unit of vorticity (in the 
usual units), because of the Z 2 divided out, and are also 
referred to as 7r-disclinations. For the case of the ferro- 
magnetic state, 7Ti(50(3)) = Z 2 , so that two nontrivial 
vortices can annihilate one another 5 . 

The other main type of topological defect is sometimes 
called a coreless vortex or skyrmion. In these the order 
exists and varies within the order parameter manifold 
everywhere in space. The topological classification re- 
quires identifying points at infinity, as if space were a 
sphere. For trivial boundary conditions (those that al- 
low the constant ordered ground state), the topological 
defects are classified by the second homotopy group 7r 2 
of the manifold. This is Z in the polar case, and trivial 
in the ferromagnetic case. However, if we wish to classify 
the vortices that carry the nonzero vorticity in a rotat- 
ing condensate, then the boundary conditions must be 
modified to allow a nonzero vorticity on the sphere (this 
modification has of the form of the Dirac string familiar 
for a magnetic monopole). In the polar case, the pres- 
ence of net vorticity may force vortices into the ground 
state. Whether the ground state contains 7r-disclinations, 
or vortices of larger vorticity each, depends on the de- 
tailed energetics, but one should note that as the total 
vorticity is an integer N v , the number of 7r-disclinations 
(more generally, the number of vortices that carry half- 



integer vorticity) must be even. On the other hand, the 
defects above that are classified by 7T2 carry no vorticity. 
For the ferromagnetic case, there are coreless vortices or 
skyrmions which are topologically-nontrivial textures in 
the order and carry nonzero vorticity quantized in inte- 
gers. We will discuss these configurations in more detail 
below. 

We now come to the borderline case in which the inter- 
action is spin-independent (SU(3) sp i n invariant), c 2 = 0. 
This is a useful starting point for small 7 also. The 
discussion is easily generalized to the case of an n- 
component order parameter with an SU(n) invariant in- 
teraction. In such a case, the target space is described 
by a complex n-component vector of unit length, pro- 
portional to the expectation value of the boson opera- 
tor, which lies in U(n)/U(n — 1) = 5 2 " -1 . For n > 1, 
7Ti = 7T2 = 0, so this space has neither point-singular vor- 
tices nor skyrmions without vorticity. For n = 1, there 
are the well-known point-singular vortices with one unit 
of vorticity each, and for n > 1 the presence of nonzero 
vorticity on the sphere induces skyrmions with integer 
vorticity in the condensate, as we will now see. 

In more detail, the Bose condensate on a sphere with 
vorticity present is described by the expectation value of 
the field operator (-0 M (z)) as z ranges over the sphere. 
Thus the possible condensates form the space of sections 
of a complex vector bundle. Such bundles are classified 
topologically (for each number of components n > 0) 
by their first Chern class, which is simply the (integer) 
number of vortices N v that we have been using, or the 
number of flux quanta in the monopole at the center of 
the sphere 2 ^. Restricting the bosons to the LLL means 
considering only the "holomorphic" sections of the same 
bundles. The LLL mean field theory performed in Sec- 
tion simply finds such holomorphic sections of lowest 
mean-field energy. Condensates in which there are N v 
vortices at which the density (magnitude squared of the 
condensate) vanishes always exist; simply take the con- 
densate entirely in one spin component. The question 
we will pursue here is the existence of coreless vortices 
or skyrmions with non-zero vorticity, in which the con- 
densate is nonzero at all points on the sphere. These 
correspond to nonsingular configurations of the order as 
described above. 

The nicest configurations of all, which serve to illus- 
trate the most elementary skyrmions, are those in which 
the density is uniform over the sphere. [Note that for 
repulsive spin-independent interactions with SU(n) sp i n 
symmetry, the ground state makes the density as uni- 
form as possible.] It is convenient to study these in the 
LLL in terms of their components b mll , fi — 1, . . . , 77, 
and m = 0, . . . , N v . These components form a matrix 
B, with to labelling rows and \i labelling columns. The 
density is uniform if and only if 

B£t a, 1 (Al) 

where I is the identity matrix (that is, is proportional 
to an isometry). Solutions to this condition exist only 
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when < N v < n — 1. For N v > n — 1, skyrmions 
will still occur, but the density will not be completely 
uniform. 

We now specialize to n = 3 again, and consider the 
special cases N v = 1 and N v = 2 for which uniform 
density condensates exist. For N v = 1, m = 0, 1 only, 
so is a 2 x 3 matrix. A solution of eq. IjAlfl for 
is a complex number times an isometry of C 2 into C 3 . 
Thus this means picking a two-dimensional subspace of 
the spin space C 3 . Notice that SU(2) orh acts by multi- 
plication of B on the left, while SU(3) sp i n acts by multi- 
plication (by the transpose) on the right, and the phase 
symmetry under U(l) acts on either side of B. Then the 
full space of solutions to eq. IjAljl . for fixed mean parti- 
cle number, is parametrized by SU (2) or b x SU (3) sp i n x 
U(1)/[SU(2) x U(l)], where the denominator represents 
the subgroup of SU(3) sp i n which has the same action as 
5t/(2)orb x U(l) on a particular solution B. This mani- 
fold is equivalent to SU(3)/%3. If we consider a quantum 
state in which bosons condense in one condensate in this 
family, then we can analyze it in terms of AT, L, L z , 
and SU(3) spm quantum numbers (p,q). It is easy to see 
that there can be no SU(3) sp [ n singlets for N v — 1 as 
construction of such a singlet requires the use of three 
distinct orbitals. Thus the broken-symmetry states can- 
not be averaged over SU(3) spin rotations to produce 
an 5/7(3) singlet (the closest one can get would be the 
BDC of Section llVl and we believe that state should be 
interpreted in this way). Nonetheless, we have shown the 
existence of configurations with single units of vorticity. 
We have now analyzed the space of solutions with uni- 
form density using the full SU (3) sp i n symmetry group. 
For general c 2 ^ 0, this symmetry is broken. In this 
case, the form of the lowest-energy solution depends on 
the energetics, and the manifold of lowest energy conden- 
sates (orbit of the solution under the broken symmetries) 
is a submanifold of that above with a lower dimension 
that depends on which solution is chosen. Results of this 
analysis have been given in Sec. 

Similarly, for N v — 2, there are three orbitals, and 
the matrix B is now 3x3. We see directly from 
eq. (|A1|I that a solution for B is proportional to a uni- 
tary matrix, and the manifold of solutions is therefore 
17(3). In terms of the symmetries present for c 2 = 
(i.e. using SU(3) sp i n symmetry), this manifold arises as 
SU(3) spin x SO(3) oth x U(l)/SO{3) x Z 3 . In this case 
the SO (3) in the denominator, which is the unbroken 
subgroup of the spin and orbital symmetries, has to be 
embedded into SU(3) as the group of 3 x 3 orthogonal 
matrices. The U(l) transformations cannot be removed 
using the orbital or spin symmetry groups. 

The uniform-density skyrmions for N v = 2 and c 2 =/= 
can be found by using a careful choice of basis. For the 
ferromagnetic case c 2 < 0, one expects that at each point 
on the sphere there should be a nonzero (in fact, largest 
possible) spin density, though its orientation varies over 



the sphere. If we use the basis of single-particle S z eigen- 
states in the sequence a =J,, 0, T, then one such solu- 
tion is given by B oc diag(l, — 1, 1). In the particular 
solution given, the spin density is J, at the north pole 
(z = in stereographic coordinates), f at the south pole 
[z — > oo), and the orientation elsewhere on the sphere 
can be found by the relation of spin and orbital rota- 
tions as described in the previous paragraph. In fact, 
the spin-density (S) itself wraps around the sphere, be- 
ing (0,0,-1) at the north pole, (0,0,1) at the south 
pole, and in the xy plane at the equator. This solu- 
tion is the 1=1 solution discussed in Sec. [V] Other 
solutions are obtained by multiplication of B by an el- 
ement of 5*0(3) x U(l). The space of these solutions 
forms the manifold SO(3) spin x SO(3) olh x U(l)/SO(3) ~ 
50(3) x U(l), a submanifold of the full (7(3) we had be- 
fore. The solution B = diag(0, —1, 0) we began with here 
was chosen to be invariant under the diagonal 50(3) sub- 
group generated by L + S. The other solutions, obtained 
by acting with either S0(3) O rb or S f O(3) sp i n , are invariant 

under a similar subgroup with generators L plus a (fixed) 
50(3) rotation of S. These solutions are very similar to 
the basic skyrmions for the spin-1/2 case, which appear 
in the v = 1 QH effect for electrons 3 ^, even though here 
for spin 1 they have N v = 2. In fact those solutions for 
a condensate of spin-1/2 bosons correspond to the n = 2 
component case with a contact interaction and SU (2) 
symmetry, as discussed briefly above. 

In the antiferromagnetic regime 7 > for N v — 2, one 
expects a different mean field ground state. As 7 — > + , 
the ground state should approach a point in a submani- 
fold of the U (3) manifold of uniform density states. Pre- 
sumably this submanifold is a different one from that for 
7 < 0, where the solutions always lie in the set we just 
described. In fact, we expect that the ground states se- 
lected at 7 — > 0+ are those with the lowest spin density 
at each point on the sphere. At large 7, the spin order 
at each point should take the polar form, with vanishing 
spin density. This presumably cannot happen with a uni- 
form density. One expects the preceding uniform density 
solutions to persist at small positive 7, and then (possi- 
bly, above a nonzero critical value of 7) for the density 
to become non- uniform. At large enough 7, we expect 
the N v = 2 solution to contain two polar vortices, each 
with vorticity one. These correspond to solutions found 
in Sec El 

For each of the manifolds that describe lowest-energy 
mean field solutions, we can consider the configuration 
space of a point moving on this manifold, and then quan- 
tize this motion (this is known as semiclassical quantiza- 
tion of the collective coordinates) . This should reproduce 
the full space of states for the N — ► 00 limit for these 
cases N v = 1, 2, but we will not go into these details 
here. 
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